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Abstract 

The  purpose  of  this  study  was  to  investigate  and  describe  die  differences 
between  a  one-degree-of-freedom  (1-DOF)  controller,  in  which  a  controller  is  placed  in 
a  single  location  in  the  control  loop,  and  a  two-degree-of>freedom  (2-DOF)  controller, 
in  vsdiich  two  separate  controllers  are  designed  ->  one  in  die  feedback  loop  and  the 
other  as  a  prefilter  in  the  forward  loop.  The  ^proach  involved  summarizing  the  major 
rules  governing  loop  shaping  for  performance  and  robustness  in  die  1-DOF  case  and 
then  extending  the  concepts  to  die  2-DOF  case.  H2  optimization  was  utilized  to 
perform  several  types  of  SISO  examples  and  one  MEMO  example  to  verify  the 
concepts.  In  all  cases,  the  2-DOF  controller  exhibited  superior  tracking  performance 
over  a  wide  range  of  plant  disturbances  and  measurement  noises  compared  to  die  1- 
DOF  model. 
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OPTIMAL  CONTROL  DESIGN  ADVANTAGES 
UTILIZING  TWO-DEGREE-OF-FREEDOM 

CONTROLLERS 

I.  Introduction 

1.1  Background 

Current  linear  control  system  design  depends  primarily  upon  the  principles  of 
either  cascade  or  feedback  compensation.  Employing  these  methods  allows  the  use  of 
such  design  techniques  as  root  locus  and  other  classical  design  techniques.  In  any 
case,  the  basic  block  diagram  used  in  the  madiematical  model  requires  that  the 
compensator  be  placed  in  a  single  location  in  die  control  loop,  as  shown  in  Figure  1- 
1(a).  Hence,  systems  diat  are  modeled  by  such  a  block  diagram  are  referred  to  as  one- 
degree^if-freedom  (1-DOF)  systems.  The  basic  control  design  objective  examined  in 
diis  thesis  is  achieving  favorable  command  following  performance  in  the  face  of  plant 
disturbances  and  sensor  or  measurement  noise.  In  1-DOF  systems,  the  single 
compensator  is  designed  to  accomplish  all  of  die  tasks  required  for  this  objective.  The 
problem  typically  encountered  with  this  type  of  system  is  that  there  is  an  inherent 
tradeoff  between  tracking  low  frequency  commands  and  rejecting  high  frequency 
measurement  noise. 

A  more  successful  w^  to  deal  with  this  problem  is  to  use  a  two-degree-of- 
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freedom  (2-DOF)  ^stem,  ^ich  is  shown  in  Figure  l-l(b).  Using  this  type  of  setup, 
dtere  are  now  two  separate  compensators  available  to  accomplish  the  design 
objectives.  The  idea  behind  this  type  of  system  is  that  the  compensator  in  die 
feedback  loop,  Ky,  takes  care  of  the  disturbance  and  measurement  noise  rejection 
objectives,  vdiile  the  prefilter,  K2,  takes  care  of  command  following  properties.  From 
die  outset,  it  is  expected  that  2-DOF  controllers  will  provide  better  overall  tracking 
performance  than  I -DOF  controllers;  two  of  the  questions  that  this  thesis  proposes  to 
answer  are  how  much  better  and  under  v^at  conditions  they  are  better. 


Figure  1-1.  (a)  1-DOF  Control  System  Block  Diagram 
(b)  2-DOF  Control  System  Block  Diagram 


1.2  Objectives 

The  primary  objective  of  diis  thesis  is  to  investigate  and  quantify  the 
differences  between  the  1<D0F  and  2-DOF  design  methodologies.  This  objective  is 
accomplished  by  observing  dre  trends  in  the  closed  loop  frequency  and  time  responses 
of  both  systems  as  plant  disturbances  and  measurement  noises  are  varied. 

1.3  Scope 

The  iq)proach  followed  is  to  first  summarize  the  major  rules  governing  loop 
shaping  for  performance  and  robusmess  in  the  1-DOF  case  and  then  extend  these 
concepts  to  die  2-DOF  case.  Then,  several  types  of  examples  are  examined  in  order  to 
illustrate  the  differences  between  die  two  systems.  The  first  set  of  examples  are 
single-innut-single-output  (SISO):  a  stable,  minimum  phase  system;  an  unstable, 
minimum  phase  system;  a  stable,  nonminimum  phase  system;  and  an  unstable, 
nonminimum  phase  system.  This  first  set  of  SISO  examples  is  explored  in  detail  to 
show  die  basic  trends  that  can  be  expected  from  each  type  of  ^stem.  The  next  set  of 
SISO  examples  involves  designing  each  ^stem  to  be  a  particular  system  type  and 
observing  its  time  responses  to  higher  order  inputs.  The  last  SISO  example  examines 
die  effect  of  low  frequency  measurement  noise  on  the  tracking  performance  of  both  1- 
DOF  and  2-DOF  systems.  Finally,  a  multiple-input-multiple-output  (MEMO)  point 
design  for  a  missile  autopilot  system  is  analyzed  to  verify  die  concepts  for  the  MIMO 
case. 

There  are  several  limitations  regarding  diis  thesis  that  are  worth  noting  here.  It 
is  not  the  purpose  of  this  thesis  to  actually  perform  designs,  in  either  the  SISO  or 
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MIMO  case.  Repres^tative  disturbance  and  measurement  noise  sizes  are  repeated  in 
all  of  the  examples  to  ^ow  the  trends  in  the  frequency  re^nse  plots  and  time 
responses  —  not  to  provide  design  data.  Also,  all  of  the  controllers  used  in  the 
examples  in  diis  thesis  were  designed  using  Hj  optimization,  \^l4lich  has  certain 
limitations,  such  as  the  two  compensators  in  the  2>DOF  system  being  constrained  to 
share  the  same  poles.  In  keeping  with  dirust  of  the  thesis  objectives,  however,  H2 
optimization  represents  a  consistent  design  mediod  with  which  the  design 
characteristics  of  l-DOF  and  2*DOF  systems  can  be  compared. 

1.4  Related  Work 

In  section  1.1  of  this  thesis,  it  was  stated  that  2-DOF  controllers  are  expected 
to  perform  better  than  1>D0F  controllers.  This  general  conclusion  has  come  about  due 
to  the  previous  published  work  of  odiers  in  this  field.  While  dte  2-DOF  controllers  in 
diis  work  are  designed  simultaneously  using  H2  optimization,  most  of  the  work  in  this 
field  has  made  use  of  other  mediods,  such  as  optimization,  or  describes  methods  of 
designing  the  two  compensators  separately  and  parameterizing  the  class  of  all 
compensators  that  satisfy  some  criteria.  For  references  utilizing  these  techniques,  the 
reader  is  referred  to  [4,  5,  8,  9,  10].  Since  the  intent  of  thesis  is  not  to  perform  actual 
designs,  but  to  provide  insight  into  the  fundamental  differences  between  1>D0F  and  2- 
DOF  systems,  other  optimization  techniques  were  not  investigated. 

1.5  Thesis  Outline 

This  diesis  consists  of  eight  chiqiters.  Chapter  I  provides  general  background 
information  on  the  potential  benefits  of  the  2'DOF  design  method  and  outlines  the 
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scope  of  dtis  work. 

Chapter  n  covers  fundamental  design  concepts  for  the  1-DOF  case,  for  bodi 
SISO  and  MIMO  cases.  Specifically,  the  attributes  which  comprise  a  "good"  loop 
shape  with  regard  to  tracking,  disturbance  rejection,  and  measurement  noise  rejection 
are  described.  System  type  concepts  are  also  reviewed.. 

Ch{4)ter  m  mirrors  the  previous  chapter  by  extending  the  ideas  of  loop  shaping 
for  performance  and  robusmess  to  the  2>DOF  case,  emphasizing  key  differences. 

Chapter  IV  provides  an  overview  of  the  benefits  and  limitations  of  the  Hj 
optimization  design  method,  which  is  used  for  all  of  the  examples  in  this  thesis. 

The  general  H2  control  problem,  as  it  is  applied  to  die  examples  in  this  diesis, 
is  described  in  Chapter  V.  In  addition  to  showing  the  appropriate  transfer  function 
matrices  and  state  spaces  for  the  1-OOF  and  2-DOF  cas^,  the  rationale  behind  all  of 
the  weight  selections  used  in  the  examples  is  provided. 

Chapter  VI  features  SISO  examples  to  illustrate  the  concepts  of  Chapters  n  and 
in.  First,  an  exhaustive  frequency  and  time  response  plot  analysis  is  performed  for  1) 
a  stable,  minimum  phase  system;  2)  an  unstable,  minimum  phase  system;  3)  a  stable, 
nonminimum  phase  system;  and  4)  an  unstable,  nonminimum  phase  ^stem.  Secondly, 
die  unstable,  nonminimum  phase  example  is  used  in  conjunction  with  altered  H2 
tracking  weights  to  compare  the  effects  of  higher  order  inputs  on  different  system 
^es.  Finally,  die  same  unstable,  nonminimum  phase  plant  is  used  to  illustrate  the 
effects  of  low  frequency  measurement  noise  on  1-DOF  and  2-DOF  systems.  This  last 
example  represents  a  prime  case  i^ere  the  advantages  of  utilizing  a  2-DOF  controller 
should  be  iqiparent. 
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Chapter  Vn  presorts  a  MIMO  missile  autopilot  example  and  explores  the 
trends  in  loop  shape  and  tracking  performance  as  the  plant  disturbance  and 
measurement  noise  are  varied. 

Chapter  Vm  provides  the  overall  conclusions  for  this  work,  as  well  as 
recommendations  for  further  study. 
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//.  Fundamental  Concepts  for  UDOF  Controllers 


The  purpose  of  this  chapter  is  to  review  the  basic  concepts  of  feedback  control 
of  linear  systems  in  the  1-DOF  case  (both  SISO  and  MIMO)  in  order  to  provide  a 
basis  for  comparison  to  the  2-DOF  case.  Chapter  m  extends  these  basic  concepts  to 
die  2-DOF  case  and  highlights  the  key  differences.  Since  the  overall  objectives  in 
control  system  design  include  command  following,  disturbance  rejection,  and 
measurement  noise  rejection,  each  of  these  concepts  will  be  examined  in  detail  ~ 
specifically,  how  the  shape  of  die  closed  loop  ^stem  should  look  in  order  to  achieve 
these  objectives.  In  addition,  die  idea  of  system  type  and  different  types  of  inputs  will 
be  compared  between  the  two  cases.  Although  the  differences  between  1-DOF  and  2- 
DOF  systems  are  best  illustrated  dirough  die  use  of  examples,  this  background  chapter 
lays  the  groundwork  for  what  should  be  expected  in  the  examples  in  die  latter 
chapters. 

2.1  Loop  Shaping 

One  of  die  most  important  aspects  of  control  system  design  is  the  concept  of 
"good"  loop  shaping.  In  most  cases,  a  control  system  designer  has  at  least  one 
objective  in  mind;  to  build  a  set  or  multiple  point  tracker,  to  reject  plant  disturbances 
or  measurement  noises,  or  to  achieve  performance  in  die  face  of  plant  uncertainty. 
Often,  all  of  these  objectives  are  important.  The  shape  of  a  favorable  loop  for  a 
closed  loop  system  depends  upon  ^hat  performance  criteria  are  desired.  This  chapter 
is  devoted  to  detailing  vdiat  characteristics  a  1-DOF  control  system  must  exhibit  to 


satisfy  each  design  objective.  A  more  complete  treatment  of  loop  shaping  is  contained 
in  [7;Ch.  4]. 

2.1.1  SISO  Case.  Figure  2-1  shows  the  standard  1-DOF,  SISO,  unity 
feedback,  cascade  compensation  block  diagram  which  is  taught  in  most  basic  control 
^stem  fundamentals  courses. 


Figure  2-1.  I -DOF  SISO  Control  System  Block  Diagram 


where 

=  reference  command  input  signal 
e  =  error  signal 
1/  =  control  input  signal 
d  =  output  disturbance  signal 
n  =  sensor  noise  signal 
y  =  output  signal 
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Relating  die  outputs  to  die  inputs  in  closed  loop  form  yields 


Ir  4- 

GK  ) 

[i  *gk)  ‘  [ 

1  *gk}  ( 

1  +  GK) 

(2.1) 


It  is  i^ipropriate  here  to  define  two  quantities  which  will  prove  important  in  the  design 

process.  First,  die  quantity  [  ^  |  (or  [/  +  GK\^  in  the  MIMO  case)  is  known  as 

( 1  *  GK) 

die  system's  sensitivity.  Second,  the  quantity  j — — J  (or  [/  +  GKY^GK  in  the 

(1  +  GK ) 

MIMO  case)  is  known  as  the  system's  complimentary  sensitivity.  Denoting  the 
sensitivity  as  S  and  the  complimentary  sensitivity  as  T,  equation  (2.1)  can  be  rewritten 
as 

y  =  Tr  *  Sd  -Tn  (2.2) 

From  this  equation,  die  iqipropriate  size  of  die  magnitude  of  GK  (denoted  \GK\  )  can 
be  determined  depending  on  the  particular  objective  desired. 

2.1.1.1  Command  Following.  For  good  tracking,  consider  first  only 
the  input  and  the  output  y.  The  transfer  function  from  to  y  is  simply  the 
complimentary  sensitivity.  It  is  ^parent  diat  if  |1  +  GK\  (and  hence  \GK\  )  is  made 
"large"  over  the  frequency  range  \^ere  commands  are  expected,  usually  low 
frequencies,  then  the  transfer  function  will  be  close  to  unity  over  those  frequencies, 
and  the  system  will  exhibit  good  tracking  characteristics. 

2. 1.1.2  Disturbance  Rejection.  Good  disturbance  rejection  is  attained 
by  examining  the  transfer  function  from  d  to  y,  which  is  equal  to  the  sensitivity. 

Given  this,  it  is  clear  that  when  (GA^j  is  large  over  the  frequency  where  disturbances 
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are  expected  (usually  low  frequencies),  the  transfer  function  become  very  small,  and 
the  disturbance  is  effectively  rejected.  Of  note  is  the  fact  that  the  requirement  that 
\GK\  be  large  tqpplies  to  both  the  command  following  and  disturbance  rejection 
performance  objectives. 

2. 1.1.3  Sensor  or  Measurement  Noise  Rejection.  Finally,  die  transfer 
function  from  n  to  y  can  be  examined  in  order  to  achieve  sensor  noise  rejection.  From 
equation  (2.1),  this  transfer  function  is  equal  to  the  negate  of  the  complimentary 
sensitivity.  In  order  to  reject  measurement  noises,  diis  quantity  must  be  made  "small" 
at  frequencies  where  noises  are  expected.  It  is  important  to  note  that  the  output  y  is 
related  to  die  command  input  and  the  sensor  noise  input  signal  n  by  the  same 
transfer  function,  the  complimentary  sensitivity.  In  the  command  following  case,  \GK\ 
needs  to  be  made  large,  and  in  the  sensor  noise  rejection  case,  \GK\  needs  to  be  made 
small,  presenting  what  could  be  a  serious  conflict  in  design  objectives.  If  the 
measurement  noise  has  any  low  frequency  content,  it  will  be  passed  dirough  to  the 
output  since  \GK\  is  large  at  low  frequencies.  Measurement  noise,  however,  is 
typically  a  high  frequency  phenomenon,  so  as  long  as  there  is  adequate  separation 
between  the  frequency  range  i^ere  most  command  inputs  fall  and  the  range  where 
measurement  noises  fall,  there  is  no  problem  in  meeting  both  design  objectives. 

2.1.1.4  Bode  Magnitude  Plot  Analysis.  The  aforementioned 
requirements  can  be  combined  together  and  represented  on  a  Bode  magnitude  plot  that 
incorporates  all  of  the  objectives  by  presenting  "barriers"  that  ]GK]  must  avoid.  Figure 
2-2  represents  such  a  Bode  plot.  The  control  system  designer  then  tries  to  find  a 
compensator  that  produces  a  stable  closed  loop  system  whose  loop  transfer  function 
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High  frequency  barrier  for 
sensor  ndse  rejectkm 


Figure  2-2.  Bode  Magnitude  Plot  of  a  "Good”  Loop  Transfer  Function 
for  a  1-DOF,  SISO  System 

|G^  satisfies  the  properties  shown  in  Figure  2-2.  In  addition  to  a  Bode  magnitude 
plot  of  the  loop  transfer  function  \GK\,  it  is  also  beneficial  to  examine  magnitude  plots 
of  the  sensitivity  function  liS]  and  the  complimentary  sensitivity  function  \7\.  These 
Bode  plots  can  be  judged  "good"  or  "bad"  according  to  the  same  requirements  that 
produce  good  command  following,  disturbance  rejection,  and  measurement  noise 
rejection. 


In  the  case  of  the  sensitivity  function,  diere  are  two  barriers  that  |5|  must  avoid 
in  order  to  meet  all  of  the  objectives.  First,  since  \GK\  has  high  gain  at  low  frequency, 
|iS]  must  be  sufficiently  small  at  low  frequency  in  order  to  achieve  a  small  steady  state 


error  in  tfie  step  response  (o&er  types  of  inputs  will  be  discussed  later  in  the  chapter). 
The  sensitivity  function  must  also  cross  the  0  dB  line  at  a  high  enough  frequency  diat 
die  system  will  have  a  sufficiently  fast  transient  response  (the  speed  of  the  system's 
transient  response  is  proportional  to  how  hi^  the  crossover  frequency  of  \GK\  is). 
Secondly,  |5|  must  not  exhibit  too  much  overshoot  above  the  0  dB  line  because  then 
the  stability  margins  begin  to  deteriorate  (stability  margins  will  be  discussed  more 
fully  later  in  this  chapter).  Combining  these  ideas  yields  a  Bode  magnitude  plot  like 
die  one  shown  in  Figure  2-3. 

A  Bode  plot  of  the  complimentary  sensitivity  function  exhibits  the  same 
information  contained  in  die  sensitivity  plot,  but  the  information  is  not  quite  as 
explicit.  Once  again,  if  |GA!|  is  large  at  low  frequency  and  small  at  high  frequency, 
dien  the  system  exhibits  good  tracking  characteristics,  disturbance  rejection,  and  noise 


Figure  2-3.  Bode  Magnitude  Plot  of  a  "Good"  Sensitivity  Function 
for  a  1-DOF,  SISO  Transfer  Function 
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rejection.  On  a  magnitude  plot  of  |7].  the  function  is  unity  (0  dB)  at  low  frequency 
and  "rolls  off  at  high  frequency  in  order  to  achieve  diese  characteristics.  The  real 
utility  of  this  plot  is  in  noting  the  point  at  wdiich  the  frequency  passes  the  -3  dB  point, 
which  is  known  as  the  system's  bandwidth,  udiich  will  be  denoted  0^.  Hie  bandwiddi 
is  a  measure  of  die  frequency  range  over  i^ich  die  system  will  track  command  inputs. 
Incidentally,  if  any  measurement  noises  fall  below  the  system  bandwidth,  then  they 
will  be  passed  dirough  to  the  output  y.  Additionally,  like  the  sensitivity  function,  the 
complimentary  sensitivity  function  must  not  exhibit  too  much  oversh'  or  the  stability 
margins  begin  to  deteriorate.  Thus,  there  is  a  robustness,  as  well  as  a  bandwidth, 
barrier  on  a  Bode  plot  of  |7].  Figure  2-4  shows  a  Bode  plot  of  a  "good"  ^stem 
complimentary  sensitivity  function. 


Figure  2-4.  Bode  Magnitude  Plot  of  a  "Good"  Complimentary  Sensitivity  Function 
for  a  1-DOF,  SISO  System 
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2.1.2  MIMO  Case.  Any  discussion  of  system  loop  shaping  would  be 
incomplete  widiout  considering  the  multivah^le  case.  Problems  of  this  type  employ 
the  same  basic  principles  as  in  the  SISO  case,  but  require  some  minor  modifications. 
Figure  2-1  is  still  applicable  to  this  discussion,  with  the  primary  difference  being  that 
all  signals  are  now  vector  signals,  and  the  plant  is  a  transfer  function  matrix  relating 
several  inputs  to  several  outputs. 

2.1.2.1  Singutar  Values.  The  first  principle  that  applies  to 
multivariable  loop  shaping  is  the  concept  of  singular  values.  In  the  SISO  case,  the 
"size"  of  a  transfer  function  is  simply  its  m^jiitude  at  the  particular  frequency  of 
interest.  When  the  magnitudes  are  plotted  over  a  range  of  frequencies,  a  Bode  plot 
results.  In  the  MIMO  case,  however,  the  plant  is  described  by  a  transfer  function 
matrix,  whose  size  cannot  be  described  by  so  simple  a  notion.  Instead,  it  is  described 
in  terms  of  its  singular  values.  The  singular  values  of  a  given  matrix,  A,  are  given  by 

vdiere 

Oj  {A)  -  singular  values  of  A 
A*  -  complex  conjugate  transpose  of  A 
(A* A)  =  eigenvalues  of  A* A 

The  minimum  singular  value  of  A,  denoted  g(i4),  is  a  measure  of  the  maximum 
attenuation  of  a  signal  by  the  transfer  function  matrix.  The  maximum  singular  value, 
o(i4),  is  a  measure  of  the  maximum  amplification  of  the  signal.  Hence,  if  a(A)  is 
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large,  tfien  die  matrix  A  is  said  to  be  large;  if  a  (A)  is  small,  then  A  is  said  to  be 

small.  As  frequency  changes,  each  individual  transfer  function  in  die  transfer  matrix 
also  changes,  just  as  in  the  SISO  case.  Therefore,  the  minimum  and  maximum 
singular  values  of  die  transfer  function  matrix  are  also  functions  of  frequency  and  can 
be  plotted.  The  resulting  plot  is  called  a  singular  value  plot  and  is  the  MIMO 
equivalent  of  a  Bode  plot. 

2.1.2.2  MIMO  Performance  Measures.  The  concept  of  singular  values 
can  be  used  to  extend  the  SISO  ideas  of  command  following,  disturbance  rejection, 
and  measurement  noise  rejection  to  the  MIMO  case.  First  of  all,  the  MIMO  equation 
relating  the  outputs  to  the  inputs  is 

y  =  [/  +  GKY^GKr^  +  [/  +  GK]'^  d  -  [/  +  GK\'^GKn  (2  4) 

Beginning  widi  an  analysis  of  command  following  performance,  the  transfer  function 
from  to  y  needs  to  tqiproach  unity  over  the  low  frequency  range  i^ere  commands 
are  expected.  This  means  that  for  all  frequencies  below  the  bandwidth  frequency, 

[7  +  GK\-^GK  «  /  V  ©  <  ©„  (2.5) 

For  disturbance  rejection,  die  sensitivity  function  must  be  minimized  at  low  frequency, 
just  as  in  die  SISO  case,  except  now  the  condition  becomes 

[7  +  G/Q-'  «  0  V  ©  <  ©0  (2.6) 

Saying  that  the  sensitivity  function  should  be  "small"  at  low  frequency,  is  equivalent 
to  saying 
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o[{/+Gi:}-'l  "sma/r 


(2.7) 


Using  the  singular  value  identity 


V  <0  <  oOp 


condition  (2.7)  becomes 


o(^) 


1 

oiA-') 


gj/  +  GK]  -large”  V  (o  <  (o^ 
vdtich,  in  turn,  may  be  simplified  to  read 


(2.8) 


(2.9) 


gJGA-]  -large-  V  o  <  ©„  (2.10) 

Once  again,  the  measurement  noise-to-output  transfer  function,  n  to  y,  is  the  negate  of 
the  complimentary  sensitivity,  so  for  all  frequencies  above  die  bandwidth  frequency 
(where  hopefully  all  of  the  measurement  noises  are). 


(/  +  G^-'GA^  «  1  V  o  >©„  (2.11) 

It  can  be  shown  that  at  high  frequency,  equation  (2.11)  holds  if 

o[G^]  -small-  V  ®  >  ©„  (2  ’2) 

Putting  all  of  this  information  together,  it  follows  that  a  singular  value  plot  of  a 
"good"  MIMO  loop  shiqie  (a[G#n)  can  be  constructed  by  satisfying  equation  (2.10)  at 
low  frequency  and  equation  (2.12)  at  high  frequency,  as  in  Figure  2-S. 

2.1.2.3  Stability  Margins,  Thus  far,  little  mention  has  been  made 
regarding  a  system's  robustness,  or  ability  to  maintain  closed  loop  stability  in  the  face 
of  plant  uncertain^^  gain  changes,  phase  changes,  or  simultaneous  gain  and  phase 
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Figure  2-5.  Singular  Value  Plot  of  a  "Good"  Loop  Transfer  Function 
for  a  1-DOF,  MIMO  System 


changes.  To  maintain  consistency,  this  thesis  will  only  be  concerned  with 
simultaneous  gain  and  phase  changes,  which  are  measured  by  independent  gain  and 
phase  margins.  This  {4)proach  could  be  considered  some\^diat  conservative,  in  that  the 
margins  will  almost  always  come  out  to  be  less  than  margins  obtained  classically  or 
one-loop-at-a-time.  However,  in  cases  w4iere  the  plant  has  imstable  poles  or 
nonminimum  phase  zeros,  simultaneous  gain  and  phase  changes  can  sometimes  easily 
destabilize  the  system,  vdiile  gain  or  phase  changes  alone  could  not.  The  independent 
gain  margin  (IGM)  and  independent  phase  margin  (IPM)  of  a  ^stem  are  given  by 
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IGM 


U  [1  -  Oo  ,  1  +  aj 


1  1 
1  +  a„  ’  1  -  tto 


/  •> 

* 

* 

IPM  = 

±2  sin"' 

“o 

.2  . 

u 

±2  sin*' 

"o 

.2. 

where 

a„=  ^S] 

flo'  "^"21/ *^-'1  =  a{[/+L]-‘L}  5in 


(2.13) 


(2.14) 


Note  that  the  variable  L  represents  die  loop  transfer  function  {GK  in  the  1-DOF  case). 
From  these  definitions  it  is  t4}parent  that  the  shape  of  the  sensitivity  and 
complimentary  sensitivity  functions  plays  a  direct  role  in  die  determination  of  stability 
margins,  in  bodi  the  SISO  and  MIMO  cases.  From  Figures  2-3  and  2-4,  the  extent  to 
i^ich  |iS|  and  |71  "bubble  up"  above  die  0  dB  line  indicates  the  amount  that  the 
margins  deteriorate.  Thus,  the  robustness  barriers  in  these  two  figures  ensure  that  die 
closed  loop  system  has  reasonable  stability  margins.  Although  no  singular  value  plot 
is  shown  for  the  MIMO  sensitivity  or  complimentary  sensitivity  functions,  die  line 
representing  the  maximum  singular  values  will  usually  exhibit  the  same  sort  of  "hump" 
above  the  0  dB  line,  and  equation  (2.13)  can  be  used  to  calculate  the  independent 
margins. 
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2.2  System  Type 

In  addition  to  defining  characteristics  of  a  linear  system's  loop  shape,  one  may 
also  be  interested  in  quantifying  its  tracking  performance  to  a  variety  of  inputs.  One 
approach  is  to  q>proximate  the  inputs  as  linear  polynomial  functions  of  time.  The 
steady  state  error  resulting  when  each  type  of  input  is  applied  to  the  system  reveals 
important  information  about  die  system.  Specifically,  for  a  stable  system,  the  degree 
of  the  polynomial  input  for  which  the  steady  state  error  is  constant  is  referred  to  as  the 
system  type.  The  following  discussion  is  excerpted  from  [3:125-132]. 

For  the  unity  feedback,  1-DOF  setup  in  Figure  2-1,  consider  the  transfer 
function  from  r^to  e  ,  which  is  equal  to  the  sensitivity  function,  or,  in  terms  of 
Laplace  variables 

e{s)  = - 1 - r(s)  (2.15) 

where  L(s)  is  equal  to  the  loop  transfer  function  G(s)K(s).  If  the  reference  input  is 
assumed  to  be  a  polynomial  of  degree  k  and  is  represented  in  the  time  domain  as 

P'«) 

where  t<.](t)  represents  die  unit  step  function,  then  the  Laplace  transform  of  the  input  is 
written 

r(s)  =  J_  (2.17) 

'  gk*\ 

Therefore,  if  Ar  =  0,  dien  the  input  is  a  step  of  unit  amplitude;  if  k-  1,  the  input  is  a 
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ramp  with  unit  slope;  \fk  =  2,  the  input  is  a  parabola  widi  unit  second  derivative;  and 
so  on. 


From  die  final  value  dieorem,  the  system  steady  state  error  is  then  defined  as 


lim  1 

-*0  ^  [1  *L(s)]s^*' 


lim  1 

fl  +Z,(^)l5* 


(2.18) 


If  the  system  is  type  0,  then  L  has  no  poles  at  die  origin  and  setting  k  equal  to  zero 
(step  input)  produces  the  following  error: 


1  +1(0)  1 


(2.19) 


vdiere  is  the  position  error  constant,  which  is  equal  to  die  DC  gain  of  L(s)  (for 
unity  feedback).  In  other  words,  if  die  loop  transfer  function  L  has  no  poles  at  the 
origin,  then  diis  error  cannot  be  zero.  Likewise,  if  the  system  is  type  I,  then  L  has  one 
pole  at  the  origin,  and  if  the  input  is  a  ramp  (k  =  1),  dien  die  steady  state  error  is 


^  lim  1  ^  lim  1  ^  1 

”  ^->0  [1  +  L{s)]s  ’  ^^0  sLis)  ” 


(2.20) 


i^ere  is  the  velocity  constant.  If  the  system  is  type  n,  L  has  two  poles  the 
origin,  and  if  the  input  is  a  parabola  (k  =  2),  dien  the  steady  state  error  is 


lim  1  ^  lim  1  ^1 

"  [1  +  L{s)]s^  s^L{s)  ~  Y^ 


(2.21) 


vdiere  Kg  is  the  acceleration  constant. 

It  is  important  to  note  that  although  a  system  of  type  k  produces  a  constant 
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error  for  inputs  whose  Laplace  transfonn  have  a  degree  of  (A  +  1),  it  perfectly  tracks 
inputs  which  have  a  degree  of  k  or  lower.  If  the  degree  of  the  input  is  {k  +  2)  or 
higher,  however,  die  system  will  produce  an  infinite  steady  state  error.  It  is  also 
important  to  note  that  the  above  results  are  only  valid  for  the  unity  feedback  case.  For 
the  non-unity  feedback  case,  the  error  constants  must  be  derived  by  solving  for  the 
reference  command-to-error  transfer  function  from  the  particular  block  diagram  and 
^plying  die  fmal  value  dieorem  as  appropriate. 
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IIL  Fundamental  Concepts  for  2-DOF  Controllers 


At  first  it  may  not  seem  necessary  to  provide  a  separate  ch^ter  on  the  basic 
linear  system  dieory  behind  2>DOF  controllers,  since  the  principles  from  which  it  is 
derived  were  just  described  in  Chapter  II.  However,  in  keeping  with  the  main  thrust 
of  this  thesis,  the  linear  system  fundamentals  as  they  apply  to  2-DOF  systems  must  be 
described  in  detail  so  that  the  expected  performance  advantages  have  a  sound  basis  in 
dieory  as  well  as  in  practice.  This  chapter  mirrors  all  of  the  concepts  covered  in 
Chapter  n,  except  the  ideas  are  extended  to  the  2-DOF  case  for  comparison  and 
contrast. 

3.1  Loop  Shaping 

In  die  2-DOF  case,  many  of  the  same  principles  of  loop  taping  from  die  1- 
DOF  case  still  apply,  but  the  added  flexibility  of  having  two  different  compensators 
manifests  itself  in  several  wiQrs. 

3.1.1  SISO  Case.  The  basic  block  diagram  for  the  2-DOF  system  is  shown  in 
Figure  3-1.  All  of  the  variables  are  the  same  as  those  defined  in  the  l-DOF  setup  in 
Figure  2-1  except  in  this  figure,  the  measured  output,  y„,  is  shown.  Of  special  note  is 
the  fact  that  the  2'DOF  block  diagram  is  derived  from  a  block  diagram  manipulation 
of  the  l-DOF  diagram,  and  the  two  are  equivalent  if  and  K2  are  equal.  If  they  are 
not  equal,  however,  then  the  equations  describing  the  system  become  quite  different 
because  is  in  the  feedback  loop,  while  K2  is  a  completely  separate  compensator, 
acting  as  a  prefilter. 
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Figure  3-1.  2-DOF  SISO  Control  Sy^^em  Block  Diagram 


Relating  the  output  to  the  inputs  in  equation  form  yields 


y 


[ 

r  + 

'  ' 

1 

d  ~ 

f 

GK, 

1  +  GK, 

V  « 

1  +  GK, 

\  J  / 

or,  more  compactly. 


y  ~  SGK^  +  Sd  -  Tn 


(3.1) 


(3.2) 


It  is  immediately  iq>parent  that  the  transfer  function  from  to  y,  which  will  be 
referred  to  as  f,  is  not  equal  to  the  complimentary  sensitivity,  which  is  still  the 

transfer  function  from  n  to  y.  Also,  the  sensitivity  and  complimentary  sensitivity 
functions  depend  only  upon  K^,  not  implications  of  these  observations  will 

become  clear  as  each  of  the  design  objectives  is  explored  in  detail. 
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3.1.L1  Disturbance  Rejection.  To  simplify  die  analysis,  the  objectives 
dealing  with  the  Ki  compensator  will  be  discussed  first.  As  in  the  1-DOF  case,  good 
disturbance  rejection  is  attained  by  examining  die  ^  to  y  transfer  function,  which  is  the 
sensitivity  function.  In  order  to  reject  a  disturbance,  die  magnitude  of  GK^,  denoted 
\GK^\,  must  be  sufficiently  large  at  low  frequency,  i^iiere  disturbances  are  expected. 

3.1.1.2  Sensor  or  Measurement  Noise  Rejection.  Measurement  noise 
rejection  is  dependent  upon  the  n  to  y  transfer  function,  which  is  equal  to  the  negate  of 
the  complimentary  sensitivity  function,  as  in  the  SISO  case.  In  order  to  reject 
measurement  noises,  \GK^\  must  be  sufficiendy  small  at  high  frequency,  vdiere  noises 
are  expected. 

3.1.1.3  Command  Following.  Thus  far,  it  has  been  shown  that  in 
order  to  attain  good  disturbance  rejection  and  measurement  noise  rejection,  restrictions 
are  placed  on  die  size  of  IGilTil  w^ich  are  ^sentially  the  same  as  those  for  the  1-DOF 
case.  However,  favorable  tracking  of  commands  is  attained  through  a  combination  of 

requirements  on  both  |GAi|  and  \GKji,  as  evidenced  by  the  f  transfer  function,  wdiich 

is  the  product  of  the  sensitivity  function,  S,  and  GK2.  This  transfer  function  must 
approach  unity  over  the  frequency  range  where  commands  are  expected  (low 
frequencies).  It  is  this  characteristic  of  the  2-DOF  model  that  sets  it  ^art  from  the  1- 
DOF  model. 

When  both  \GKi\  and  nearly  equal  to  each  odier,  the 

system  will  track  commands  that  are  in  the  frequency  range  of  interest.  However,  it  is 
not  necessary  for  \GK^  to  be  large  at  low  frequency  in  order  to  achieve  good  tracking 
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—  it  only  needs  to  be  large  enough  to  reject  the  specific  disturbance  expected.  |Cr^2l 
can  be  adjusted  to  have  its  magnitude  be  whatever  it  needs  to  be  to  make  f  equal  to 
unity. 

3.1.L4  Bode  Magnitude  Plot  Analysis.  The  requirements  just 
mentioned  present  some  interesting  possibilities  in  the  way  of  Bode  plots.  Aldiough 
there  are  countless  combinations  of  different  command,  disturbance,  and  measurement 
noise  sizes,  the  two  cases  presented  here  are  low  frequency  commands  with:  1) 
moderate-sized  high  frequenQr  noise  with  a  large  low  frequency  disturbance;  and  2) 
moderate-sized  high  frequency  noise  with  a  small  low  frequency  disturbance. 

The  plant  in  all  cases  is  assumed  to  be  stable  and  minimum  phase. 

In  the  first  case,  the  large  disturbance  is  the  driving  factor  in  determining  the 
shape  of  the  loop  transfer  functions.  At  low  frequency,  \GKi\  must  be  large  in  order  to 
reject  the  disturbance,  ^^le  at  high  frequency,  die  measurement  noise  is  large  enough 
to  drive  down  die  size  of  \GK^.  Meanw^iile,  \GKji  matches  |GAr,|  almost  identically  at 
low  frequency,  but  vdien  starts  to  roll  off  afrer  it  passes  the  bandwidth 
frequency,  \GK<^  is  able  to  level  out  at  the  0  dB  line  in  order  to  maintain  near-perfect 
tracking  out  to  a  much  higher  frequency.  The  trends  just  described  for  this  first  case 
are  shown  in  Figure  3-2. 

In  the  second  case,  the  disturbance  is  not  so  large  anymore,  but  the  size  of  the 
measurement  noise  is  still  moderately  large,  so  the  measurement  noise  exerts  the  most 
influence  on  die  shape  of  the  loop  transfer  functions.  The  overall  effect  of  a  dominant 
measurement  noise  is  to  hold  down  \GK^\  as  much  as  possible  over  all  frequencies. 


Low  frequency  barrier  for  good 
disturbance  rejection  IGK|I 


High  frequency  barrier  for 
sensor  noise  rejection  IGK,I 


Figure  3-2.  Bode  Magnitude  Plot  of  "Good”  Loop  Transfer  Functions 
for  a  2-DOF,  SISO  System  (Large  Disturbance  Case) 

especially  high  frequencies.  will  still  be  sufficiently  large  at  low  frequency  to 
reject  the  disturbance,  however.  accommodates  diis  situation  by  remaining  at  or 

around  0  dB  over  as  wide  a  frequency  range  as  possible,  diereby  keeping  f  equal  to 

unity  and  achieving  good  tracking.  Figure  3-3  shows  the  trends  just  described  for  die 
second  case. 

Another  interesting  point  regarding  the  2-DOF  setup  is  the  shape  of  the 
sensitivity  function.  Recall  from  Figure  2-3  that  in  the  1-DOF  case  there  were  two 
separate  barriers  on  the  sensitivity  function  Bode  plot  which  corresponded  to  a  steady 
state  error  and  speed  of  response  boundary  and  a  robustness  boundary.  As  previously 
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Figure  3-3.  Bode  Magnitude  Plot  of  "Good"  Loop  Transfer  Functions 
for  a  2-DOF,  SISO  System  (Small  Disturbance  Case) 


mentioned,  diere  is  no  low  frequency  barrier  which  corre^nds  to  good  command 
following  in  a  2-DOF  system.  Instead,  |G^||  must  only  be  large  enough  to  reject  die 
l^lpropriate  size  disturbance  encountered.  In  fact,  a  closer  examination  of  the  2-DOF 


sensitivity  term. 


1 

1  +  GK, 


,  reveals  that  if  the  quantity  GK^  takes  on  a  very  small  (or 


negative)  value,  then  liSI  could  take  on  valu^  ranging  between  0  and  oo  (or  -oo  to  qo 
dB).  Most  practical  or  realizable  plants  and  control  systems  would  never  utilize  many 
of  these  |5|  values  due  to  poor  stability  margins,  but  if  this  were  not  a  consideration, 
then  those  values  would  dieoretically  be  possible.  Therefore,  the  only  consideration 
remaining  for  the  2-DOF  sensitivity  function  Bode  plot  is  the  robustness  barrier  on  the 
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maximum  magnitude  of  |iS|.  The  2>D0F  |iS|  plot  resembles  the  l>DOF  |iS|  plot  in  diat 
die  curve  evratually  steadies  out  at  the  0  dB  line,  but  it  does  not  necessarily  approach 
the  0  dB  line  from  low  frequency.  Figure  3-4  shows  the  results  of  the  foregoing 
analysis  of  die  2-DOF  sensitivity  function. 

In  a  similar  fashion,  an  examination  of  die  2-DOF  complimentary  sensitivity 


function. 


GK. 


1  +  GK, 


,  reveals  that  if  the  quantity  GK^  takes  on  very  small  (or 


negative)  values,  then  |7|  could  also  take  on  values  ranging  between  0  and  qo  (or  -oo 
and  00  dB).  Once  again,  deteriorating  stability  margins  are  the  main  consideration  in 
the  shape  of  the  2-DOF  complimentary  sensitivity  function  Bode  magnitude  plot.  It 
should  be  noted  diat  even  in  die  2-DOF  setup,  disturbance  rejection  and  measurement 


Figure  3-4.  Bode  Magnitude  Plot  of  Possible  "Good"  Sensitivity  Functions 
for  a  2-DOF,  SISO  System 
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noise  rejection  are  competing  objectives.  The  plot  will  still  maintain  Ae  general  shape 
of  the  1-DOF  complimentary  sensitivity  plot  and  roll  off  at  some  bandwidth  cutoff 
frequenQT,  but  the  overall  magnitude  at  low  frequency  is  quite  variable  in  the  2-DOF 
case.  The  resulting  2-DOF  complimentary  sensitivity  function  from  this  analysis  is 
shown  in  Figure  3-S. 

Anodier  idea  to  reiterate  here  is  die  fact  that  the  complimentary  sensitivity 
function  does  not  provide  information  about  the  frequency  range  over  which  a  2-DOF 
system  will  track  commands.  The  2-DOF  equivalent  to  the  1-DOF  input-to-ouq>ut 
transfer  function  is  SGK2.  This  transfer  function  needs  to  be  equal  to  unity  over  as 
wide  a  frequency  range  as  possible  to  provide  favorable  tracking.  Since  GK^  and  GK2 


Figure  3-5.  Bode  Magnitude  Plot  of  Possible  "Good"  Complimentary  Sensitivity 
Functions  for  a  2-DOF,  SISO  System 
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are  independent,  a  decrease  in  the  size  of  one  can  be  offset  by  an  increase  in  die  size 
of  die  odier  in  order  to  maintain  tracking  performance.  This  is  die  essence  of  having 
two  degrees  of  freedom.  Therefore,  a  "good"  Bode  magnitude  plot  of  \SGKji  is 
identical  to  die  one  shown  for  the  1-DOF  case  in  Figure  2-4. 

3.1.2  MIMO  Case.  The  extension  of  die  2-DOF  SISO  loop  shaping  ideas  to 
the  MIMO  case  is  relatively  straightforward,  making  use  of  singular  values  to  define 
the  sizes  of  transfer  function  matrices.  Starting  with  the  overall  r^to  y  transfer 
function,  equation  (3.1)  becomes 

>»=[/+  GKX^GK^r^  +  [/  +  GKX^d  -  [/  +  GKy^GK^n  (33) 

Equation  (3.3)  can  also  be  written  in  terms  of  T,  and  f ,  as  in  equation  (3.2).  Good 

disturbance  rejection  requires  that  die  sensitivity  function  be  sufficiently  small  at  low 
frequency  to  accommodate  any  disturbances  expected,  v^ich  can  be  expressed  as 

g  \GK^  '^sufficiently  large*  V  ©  <  ©q  (3.4) 

Conversely,  good  measurement  noise  rejection  requires  that  |7]  be  small  at  high 
frequency,  v^ich  can  be  written  as 

a[GJir,]  "sufficiently  smalF  V  ©  >  ©^  (3-5) 

As  far  as  MIMO  command  following  performance  goes,  once  again  the  output 
is  dependent  upon  the  behavior  of  both  Ki  and  Kj.  In  the  case  where  the  disturbance 


is  relatively  large,  tiien  g[G^,]  will  be  large  at  low  ftequracy  in  order  to  reject  die 
disturbance,  and  otG/sT,]  will  also  be  large  at  low  frequency  to  matcli  g[G^,].  At 
high  frequency,  however,  o[G^,]  must  roll  off  in  order  to  reject  die  qipropriate 

measurement  noise.  Tracking  is  maintained  by  gLG/Tj]  staying  at  0  dP  out  to  a 
higher  frequency.  In  the  case  where  the  disturbance  is  relatively  small  and  the 
measurement  noise  is  relatively  large,  then  o[GAr,]  may  well  be  below  0  dB  over  all 

frequencies.  In  that  case,  a[GKj\  must  remain  at  0  dB  from  low  frequency  out  to  as 

high  a  frequency  as  possible  to  provide  satisfactory  tracking  requirements. 

The  MIMO  sensitivity  and  complimaitary  sensitivity  functions  essentially 
follow  the  same  rules  as  they  do  in  the  SISO  case.  Both  singular  value  plots  have 
only  robustness  barriers,  with  no  barriers  regarding  speed  of  response,  stearfy  state 
enor,  or  bandwidth.  The  maximum  singular  values  of  both  |5I  and  |71  must  stay  below 
the  robustness  barrier  in  order  for  the  ^stem  to  maintain  adequate  stability  margins. 
There  are  no  figures  shown  for  the  2-DOF  MIMO  loop  sh^es,  but  it  is  hoped  that  the 
foregoing  discussion  provides  adequate  insight 

Overall,  there  are  several  marked  differences  between  the  loop  shaping 
characteristics  of  a  l-DOF  and  a  2-DOF  system,  the  key  one  being  in  the  tracking 
requirements.  In  the  1-DOF  system,  the  tracking  and  disturbance  rejection 
requirements  dictate  tiiat  \GK\  be  large  at  low  frequency,  while  it  must  be  small  at  high 
frequency  in  order  to  reject  any  measurement  noise.  These  two  competing 
requirements  limit  how  well  a  1-DOF  system  can  track  because  |GA)  can  only  have  a 
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large  magnitude  out  to  a  certain  frequency  before  it  must  start  to  roll  off.  In  die  2- 
DOF  case,  diere  still  exists  the  same  conflict  betweai  disturbance  and  measurement 
noise  rejection  requirements  (this  time  with  \GKi\  )  but  die  ^stem's  tracking  ability  is 
not  exclusively  tied  to  |G^i|.  \GKji  provides  the  needed  flexibility  in  die  system  by 
being  free  to  be  whatever  magnitude  it  needs  to  be  in  order  to  follow  commands  out  to 
a  much  higher  frequency.  In  fact,  the  example  presented  in  Chapter  VI  of  this  diesis 
puts  that  added  flexibility  to  the  test  by  introducing  measurement  noise  failing  in  the 
same  low  frequency  range  that  commands  and  plant  disturbances  do.  Hie  2'DOF 
system,  therefore,  can  be  considered  to  have  two  main  parts;  an  "inner  loop",  which 
takes  care  of  noise  and  robusmess  requiremoits,  and  an  "outer  loop"  (even  though  it  is 
not  really  a  loop  at  all),  u^ch  takes  care  of  tracking  requirements.  These  conclusions 
^piy  to  both  the  SISO  and  MIMO  cases,  as  appropriate. 

3,2  System  Type 

For  a  2-DOF  system,  or  any  non-unity  feedback  ^stem  for  that  matter,  the 
simple  rules  governing  system  type  and  the  computation  of  the  error  constants  in  the 
1-DOF  unity  feedback  case  do  not  apply.  Hie  basic  principle  is  die  same,  however. 
Once  again,  from  the  basic  2-DOF  block  diagram  from  Figure  3-1,  die  analysis  begins 
with  the  error  signal  r^(s)  -  y(s).  In  die  1-DOF  case,  the  transfer  function  between  this 
signal  and  the  reference  command  was  equal  to  the  sensitivity  function,  which  made 
the  subsequent  calculations  more  straightforward.  In  die  2-DOF  case,  however,  there 
is  no  signal  direcdy  on  die  block  diagram  vsdiich  represents  the  quantity  r^(s)  -  y(s),  so 
one  must  be  created.  Figure  3-6  is  a  modified  version  of  Figure  3-1  that  shows  the 
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error  signal  e^qplicitly. 


Figure  3-6.  2-DOF  Block  Diagram  Showing  Error  Signal 


From  this  modified  block  diagram,  the  error  signal  e  can  be  ^own  to  be 


e(j)  =  -  ^(s)  =  ^(0  - 


GK. 


(I  ^GK,) 


rM 


or 


(3.6) 


e(s)  =[1  -F{s)]r^{s) 


(3.7) 


where  F(s)  is  defined  as 


GK, 

F(s)  = 


1  +GKj 


(3.8) 


Since  F(sJ  contains  both  and  IC2,  it  is  clear  that  the  error-to-reference  command 
transfer  function  will  not  have  a  clean,  clear-cut  solution  like  the  1-DOF  case. 
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Instead,  die  expression  for  die  steady  state  error  becomes 


Urn  [1  -F{s)]  ^  lim  [1  -  Fjs)] 

"  s  — >0  s*  * '  ^  s  * 


(3.9) 


Once  again,  die  system  type  is  die  smallest  value  of  A  at  which  die  steady  state 
error  is  a  nonzero  constant  (not  infinity).  The  input  tracking  characteristics  for  the  2- 
DOF  system  are  also  the  same  as  the  1>D0F  case.  If  a  system  is  type  0  and  die  input 
is  a  step,  then  a  constant  steady  state  error  will  result.  If  die  system  is  type  I,  then  it 
will  track  a  step  perfecdy  and  produce  a  constant  steady  state  to  a  ramp  input.  If  the 
system  is  type  n,  dien  it  will  track  a  step  and  a  ramp  perfecdy  and  produce  a  constant 
steady  state  enor  to  a  parabolic  input.  The  same  logic  applies  to  systems  of  greater 
type  and  inputs  modeled  by  higher  order  polynomials.  The  primary  difference 
between  1-DOF  and  2-DOF  systems  widi  regard  to  system  type  is  diat  in  the  2-DOF 
case,  one  cannot  necessarily  attribute  die  number  of  poles  or  zeros  at  the  origin  in 
either  GKi  or  GK2  to  the  determination  of  the  system  type,  as  in  the  1-DOF  case. 
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IV.  H2  Optimization 


In  comparing  and  contrasting  the  single-  and  two-degree-of-freedom  design 
methodologies,  it  is  iqjpropriate  to  use  one  design  technique  that  is  consistent  between 
the  two.  optimization  represents  a  relatively  straightforward  design  technique 
\^ich,  because  it  is  an  optimization  procedure,  produces  die  "best"  compensator  for 
the  particular  model  specified.  It  is  particularly  appropriate  here,  as  the  system  is 
modeled  with  noises  acting  on  it.  It  also  allows  for  direct  comparison  between  the 
two  mediodologies  under  varying  input  parameters.  The  following  summary  of 
optimization  was  largely  taken  from  [6]  and  [2],  i^ich  should  be  consulted  for  a  more 
detailed  description  of  diis  method. 

The  objective  of  H2  optimization  is  to  produce  a  compensator  which  minimizes 
the  2-norm  of  a  transfer  function  or  transfer  function  matrix  from  exogenous  inputs 
(commands,  disturbances,  or  measurement  noises)  to  controlled  ou^uts.  The  standard 
design  setup  is  shown  in  the  diagram  in  Figure  4-1,  with  w  representing  the 
exogenous  inputs  (assumed  to  be  zero  mean,  unit  intensity  white  Gaussian  noises)  and 
z  representing  the  controlled  outputs.  P  is  die  plant,  which  includes  all  system 
dynamics  and  weights  on  the  exogenous  inputs  and  controlled  outputs.  K  is  the 
compensator,  u  is  the  controlled  input  to  die  plant,  and  y  is  the  measured  plant  outputs. 
P  can  be  partitioned  in  such  a  way  that 

z  ^  P„  ^tu  w  (4  1) 

yi  [  J  I  "  . 
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w 

u 


Figure  4>1.  H2  Design  Diagram 

The  closed-loop  transfer  function  is  die  transfer  function  from  w  to  z  and  can  be 
written  as  a  linear  fractional  transformation  (LFT)  of  P  and  K  from  Figure  4-1  as 


PMi  -  p^Krp^ 


(4.2) 


The  H2  design  procedure  involves  designing  a  stabilizing  compensator  K  that 
minimizes  the  2-norm  of  z,  vdiich  is  equivalent  to  minimizing  the  2-norm  of  T„,  given 
by 


inf 

K  adm 


inf 

K  adm 


(4.3) 


where 


±  J  /r  [r„-Oa»  da 


(4.4) 
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Since  tiie  H2  algorithms  utilized  in  this  thesis  employ  state  space  methods,  it  is 
first  necessary  to  represent  the  plant  differential  equations  in  state  space  form.  The 
state  space  representation  of  the  nine  matrices  that  comprise  the  P  matrix  are 


i  =  Ax  +  Bw  +  Bu 


z  =  C  X  +  D  w  +  Du 


y-Cx+Dw+Du 

y  yw  yu 


>Adiich  can  be  written  in  transfer  matrix  form  as 


(4.5a) 

(4.5b) 

(4.5c) 


P  = 


A 

B  B 

W  U 

D  D 

c. 

D  D 

(4.6) 


There  are  several  assumptions  regarding  these  nine  matrices  which  are  inherent 
to  the  well-posedness  and  solvability  of  the  H2  control  problem; 

(i)  =  0 

(ii)  =  0 

(iii)  {A,  B^)  stabilizable  and  (C^  A)  detectable 

(iv)  DyJ)yJ  full  rank 


(V) 


(vi) 


A  -j&I  B^ 
C  D 

t  XU 

A  -  j(Ql  B^ 


D.. 


yw 


has  full  column  rank  for  all  <0 


has  full  row  rank  for  all  go 
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Condition  (i)  arises  due  to  the  fact  that  the  transfer  function  2>norm  is  also 
proportional  to  die  area  under  the  singular  value  plot  of  vdiich  must  be  finite.  If 
is  nonzero,  then  die  closed  loop  transfer  function  will  be  infmite.  Condition  (ii), 
on  is  not  necessary  but  facilitates  an  easier  development  of  a  solution.  Condition 
(iii)  is  necessaiy  for  stabilizing  compensators  to  exist.  Condition  (iv)  is  required  so 
that  there  is  a  penalty  on  control  usage  and  so  that  there  are  no  perfect  measurements; 
otherwise,  a  singular  control  problem  resulte.  Conditions  (v)  and  (vi)  ensure  die 
existence  of  stabilizing  solutions  to  the  algebraic  Riccati  equations  that  ^pear  in  the 
problem  solution.  Condition  (iv)  can  be  strengthened  to  =  I  and  DyJDyJ  =  / 

in  order  to  simplify  the  final  form  of  the  equations  by  scaling  u  and  y.  A  detailed 
explanation  of  this  scaling  process  is  contained  in  [7:70-71].  The  H2  solution  that 
follows  assumes  that  u  and  y  have  been  properly  scaled;  scaled  quantities  will  be 
shown  with  a  tilde  (~). 

If  suboptimal  compensators  are  considered,  die  family  of  all  compensators  that 
satisfy  |7’^l2  ^  ot,  with  a  ^  Oq,  is  given  by  die  LFT  of  J(s)  and  the  constrained 
freedom  parameter  Q(s),  shown  in  Figure  4-2.  In  the  suboptimal  case,  the  constraints 
on  Q(s)  are  ^  €  RHj  and  \Q\2  ^  RHj  is  the  space  of  all  stable,  stricdy 

proper  transfer  functions  with  real,  rational  coefficients.  For  the  case  vsdiere  Q(s)  is 
equal  to  zero,  the  transfer  function  J„y(s)  is  equal  to  die  optimal  H2  compensator 
a  nonzero  Q(s)  parameter  produces  suboptimal  compensators.  The  transfer  function 
J(s)  is  given  by 
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K 


Figure  4-2.  (J,  Q)  Parameterization  Diagram  for  H2  Suboptimal  Compensator 
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The  X2  and  Yj  matrices  are  die  real,  unique,  symmetric  positive  semidefmite  solutions 
to  the  algebraic  Riccad  equations 


and 


(A  -B,d’c,)’X,  *X,(,A  -BfilC,)  -X,Bj,%  *  C?/d,  -0  (■»  >•) 


where 


and 


D„D’)C, 


6  -  B  (I  -  d’d  ) 

"tr  w\  yw 


(4.13) 

(4.14) 


The  resulting  reverse  scaled  optimal  compensator  state  space  is  given  by 


v^ere 


l£r 

< 

10 

s;'K^ 

(4.15) 


(4.16a) 

(4.16b) 


The  final  item  of  interest  regarding  Hj  optimization  centers  around  the  order  of 
the  compensators  produced  by  the  algorithms.  In  general,  the  H2  algorithm  produces  a 
compensator  that  has  the  same  order  as  the  plant,  plus  the  number  of  states  contained 
in  the  (fynamic  weights  used  to  shape  the  closed  loop  system  (weight  shapes  will  be 
discussed  in  Chiqiter  V).  It  should  be  noted  that  in  the  2-DOF  case,  the  H2  algorithm 
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produces  two  compmsators  that  share  die  same  A  matrix  and  dius,  the  same  poles  (die 
zeros  of  each  may  be  different,  however).  The  fact  that  Ky  and  must  share  die 
same  poles  could  be  considered  a  limitation  to  utilizing  die  method  of  designing  both 
compensators  simultaneously.  However,  all  of  die  examples  in  this  thesis  make  use  of 
the  H2  method  of  designing  bodi  compensators  simultaneously  because  this  mediod 
still  illustrates  the  advantages  of  using  2-DOF  controllers  when  tracking  is  one  of  the 
primary  control  objectives.  Furthermore,  under  the  assumption  that  the  exogenous 
inputs  are  \ydiite  Gaussian  noise  and  the  ou^ut  energy  is  to  be  minimized,  this 
controller  order  is  optimal. 
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K  General  Problem  Formulation 


The  puipose  of  tiiis  chapter  is  to  present  the  setup  of  the  general  H2  problem 
for  die  1-DOF  and  2-DOF  cases,  along  widi  all  of  the  assumptions  inherent  to  the 
problem  and  the  rationale  for  selecting  die  dynamic  weight  shapes  for  the  H2 
optimization  algoridim.  The  analysis  begins  with  a  state  space  description  of  each 
system,  which  provides  all  of  the  variables  necessary  for  die  P  matrix.  Since  H2 
optimization  minimizes  the  2-norm  of  T^,  or  die  energy  of  selected  outputs  to  i^te 
noise  inputs,  die  outputs  chosen  for  the  z  vector  are  not  arbitrary.  These  outputs  must 
be  representative  of  the  "errors"  to  be  minimized.  In  general,  dynamic  weights  reflect 
the  desired  frequency  range  of  an  input  that  gets  fed  through  to  the  corresponding 
output.  All  of  diese  factors  will  be  discussed  in  this  chrqiter,  and  die  dynamic  weights 
chosen  to  illustrate  design  trends  in  the  latter  examples  will  be  shown. 

5.1  1-DOF  Case 

The  1-DOF  block  diagram  used  for  die  H2  formulation  is  shown  in  Figure  S-1. 
For  this  figure,  Wj  is  the  disturbance  signal  and  w„  is  the  measurement  noise  signal. 
The  governing  control  law  for  H2  optimization  is  u  =  Ky,  with  all  of  the  transfer 
functions  relating  the  outputs  to  the  inputs  being  derived  from  this  law.  The  most 
^ticeable  difference  between  Figure  S-1  and  Figure  2-1  is  the  presence  of  weighting 
blocks  in  the  block  diagram  of  Figure  S-1.  As  mentioned  in  the  introduction  to  this 
chapter,  the  weights  help  in  the  design  of  a  compensator  vdiich  produces  a  system 
with  a  closed  loop  shape  that  suits  the  desires  of  the  control  system  designer.  It  is 


Fifure  5-1.  1-DOF  Optimization  Block  Diagram 


easy  to  see  how  the  different  weights  affect  different  aspects  of  die  system  by 
examining  die  transfer  function  matrix  between  the  outputs  (ideal  output  (zi),  control 
usage  (Z2),  and  tracking  error  (zj))  and  die  exogenous  inputs  (commands  (r^), 
disturbances  (w and  measurement  noises  (w„)).  It  is  given  by 


WJ  W^SW^  -W^TWj 

= 

W^KS  -W^KSW^  -W^KSW^ 

^3. 

W^S  -W^SW^  WJW^ 

If 

"J 

Incidentally,  minimizing  the  transfer  function  between  and  Zj  is  not  a  control 
objective,  since  it  is  desired  for  it  to  be  unity  for  favorable  tracking.  However,  Zj  is 
included  in  the  Hj  setup  so  that  the  to  z^  and  to  Zj  transfer  functions  can  be 
minimized. 

The  parameters  of  interest  in  most  control  system  designs  are  the  sensitivity,  S, 


and  the  complimentary  sensitivity,  T.  From  die  given  transfer  function  matrix,  S  is 
affected  by  IFj,  and  W^,  vdiile  T  is  affected  by  W^,  Wy  and  Wy  At  first  glance,  it 
may  seem  as  if  these  weights  would  conflict  widi  each  other  if  diey  were  set  to 
differoit  values,  and  as  it  turns  out,  this  is  true.  However,  diere  is  a  logictd  mediod 
for  choosing  the  weights  such  that  dieir  valu^  make  intuitive  sense  and  the  resulting 
compensator  produces  a  good  closed  loop  shape. 

The  state  space  representation  of  Figure  S-1  provides  some  insight  into  the 
weight  selection  rationale.  The  nine  matrices  comprising  the  P  matrix  are  given  by 


A  = 


= 


C  = 


0 

0 

0  0 

0 

A, 

0 

0  B,C 

4  0 

0  0 

A 

0  0 

0 

-B,C^  0 

0 

A,  -B,C,  0 

0  0 

0 

0  A, 

0 

0  0 

0 

0  0 

^5 

)  0 

0 

B 

g 

)  B,D, 

0 

B^D 

)  0 

0 

B  = 

2 

'a 

0 

-B,L 

>  ^4 

0 

0 

)  0 

^5, 

0 

c, 

c 

1  0  D,C,  0 

0  0 

c 

'j  0  0 

1 

0 

c 

1  C,  -DjC^  0 

-C  0  0 

g 

0 

-c  -c 

^4  '^5 

1 

(5.2a) 


(5.2b) 


(5.2c) 
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0 

0 

D  « 

IW 

0 

0 

0 

II 

,^3 

-D3D, 

0 

d„-[i-d.-d,]  Z)„-(-i).l 

The  state  vector  is  given  by 


(5.2d) 


plant  states 

controlled  output  weight  states 

control  usage  weight  states 

*3 

tracking  weight  states 

disturbance  weight  states 

.*5  . 

measurement  noise  weight  states 

If  the  assumptions  discussed  in  Chapter  IV  regarding  H2  optimization  are  applied  to 
the  above  state  space,  dien  certain  characteristics  of  some  of  the  weights  can  be 
determined,  such  as  whedier  diey  are  strictly  proper  or  just  proper.  To  begin  widi, 
must  equal  zero,  so  all  of  die  terms  in  that  matrix  must  be  set  to  zero,  including  £>3 
and  Dj  and/or  D4.  Since  is  equal  to  zero,  then  Dg  =  0,  which  means  that  die  plant 
must  be  strictly  proper.  Since  D^Jd^  must  have  full  rank,  then  Dj  must  be  nonzero 
(actually  D2D2  must  be  full  rank).  The  fact  that  DyJ)yJ  must  have  full  rank 
indicates  that  and  may  or  may  not  be  equal  to  zero.  For  all  of  die  examples 
illustrated  in  this  thesis,  the  following  assumptions  will  be  made  regarding  the 
dynamic  weights: 

(i)  Wj  is  a  constant  (Dj  nonzero,  with  no  x,  states  and  thus  no  A^,  B,,  or  C,) 
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(ii)  Wj  is  a  constant  (Dj  nonzero,  with  no  Xj  states  and  thus  no  A2,  B2,  or  C2) 

(iii)  W3  is  represented  by  a  strictly  proper  transfer  function  (D3  =  0) 

(iv)  W4  is  represented  by  a  strictly  proper  transfer  function  (£>4  =  0) 

(v)  W3  is  represwted  by  a  (nonstrictly)  proper  transfer  function  (£>;  nonzero) 
If  these  assumptions  are  i4>plied  to  the  state  space  givra  in  the  equations  of  (S.2),  die 
new  state  space  becomes 


0  0  0 
-B,C^  A,  -B,C,  0 


0 

0 


0 

0 


0 

A. 


0  0  0 

0 

0 

B  « 

w 

0  B^  0 

0  0/?, 

B  «= 

0  D,C,  0 
0  0  0  0 
0  Cj  0  0 


0  0  0 

0 

- 

0  0  0 

D  ® 

n 

p 

0  0  0 

0. 

=  [-C,  0  -C,  -C,]  0  -D,] 


D_ 


[0] 


(5.3) 


Whenever  a  weight  is  assigned  to  be  a  constant,  as  in  the  case  of  and  IPj,  only  the 
D  term  remains,  and  diat  state  can  be  removed  from  die  original  state  space.  Also,  W2 
is  a  measure  of  the  amount  of  control  power  used  in  the  ^stem,  and  it  is  given  die 
symbol  p  (a  diagonal  matrix)  to  show  that  it  can  be  "tuned”  to  speed  up  or  slow  down  the 
system. 

5,2  2-DOFCase 

The  2-DOF  block  diagram  used  in  the  optimization  algorithm  is  shown  in 
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Figure  S-2.  All  of  the  weights  are  located  in  die  same  places  as  die  1-DOF  ^stem, 
although  block  diagram  interpretation  yields  a  slighdy  different  transfer  function 
matrix  relating  die  z  outputs  to  the  exogenous  inputs.  It  is  givra  by 

W^SGK^  -WJW^]r  1 

Zj  -  W^[I  -  K^SG]K^  -W^K^SW^  (5  4) 

Wj[l  -  SGKJ  -W^SW^  WJW^  K. 

Once  again  the  sensitivity  and  complimentary  sensitivity  functions  are  affected  by  die 
same  weights  as  in  the  1-DOF  case.  Thus,  the  same  approach  to  weight  selection  will 
be  used  in  this  case.  The  primary  difference  betwera  the  1-DOF  and  2-DOF  cases  is 
the  implementation  of  the  control  law  u  =  Ky.  In  the  1-DOF  case,  diere  was  a  single 
compensator,  K,  and  only  one  control  input,  y.  For  the  2-DOF  system,  the  control  law 
must  be  modified  to  read 


Figure  5-2.  2-DOF  H2  Optimization  Block  Diagram 


u  «  Ky  ^  K^] 

As  discussed  in  Chapter  m,  die  addition  of  the  second  compensator  changes  die  input- 
output  relationships  and  accounts  for  the  differences  between  die  1-DOF  and  2-DOF 
transfer  function  matrices  in  (S.l)  and  (5.4). 

The  2-DOF  state  space  representation  of  Figure  S-2  is  actually  quite  similar  to 
die  I'DOF  state  space,  the  only  differences  being  in  die  matrices  relating  to  die 
control  input,  y  (Cy,  D^,  and  DyJ.  With  all  of  die  assumptions  regarding  Hj 
optimization  taken  into  account,  die  final  2-DOF  state  space  is  given  by 


y. 


(5.5) 


.4^  0  0  0 

0  0  0 
0  0  0 


o 

o 

o 

o 

o 

1 

0 

0  B,  0 

K  = 

0 

0  0  B, 

0 

£>,C^  0  0 

0  0  0  0 

0  Cj  0  0 


COO 
0  0  0 
0  0  0 


C,  0  -c,  -c, 

0  0  0  0 


0  0  -D, 
7  0  0 


£)_ 


0 

0 


Since  there  is  no  substantial  change  in  the  2-DOF  state  space  regarding  the  D  terms  of 
the  dynamic  weights,  all  of  the  assumptions  regarding  vdiether  a  weight  is  stricdy 
proper,  proper,  or  equal  to  a  constant  (assumptions  (i)  through  (v)  in  Section  5.1)  are 
still  tqiplicable  from  the  1-DOF  case. 
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5.3  Weight  Selection 


As  far  as  choosing  actual  values  for  the  weights,  diere  are  basically  two 
schools  of  diought.  First,  a  weight  can  be  viewed  as  a  filter  through  v^ch  only  the 
frequ«tcies  of  interest  are  allowed  to  pass.  This  logic  is  particularly  applicable  to 
noise  weights.  Since  all  of  die  disturbances  and  measurement  noises  in  H2 
optimization  are  assumed  to  be  white  Gaussian  noises,  it  makes  sense  to  only  pass 
duough  certain  frequency  ranges  of  interest  to  die  system.  Second,  a  weight  can  be 
viewed  as  a  penalty  on  a  signal  over  a  particular  frequency  range.  In  other  words,  die 
weight  is  chosen  in  such  a  way  that  its  magnitude  plot  looks  like  the  inverse  of  what 
the  designer  wants  the  signal  of  interest  to  look  like.  It  is  possible  for  these  two 
weight  selection  qiproaches  to  conflict  widi  each  odier,  but  in  general,  they  are  not 
mutually  exclusive.  The  weights  chosen  in  this  section  are  typical  of  many  aircraft 
control  problems  and  will  be  used  in  all  of  the  examples  in  diis  thesis. 

Beginning  with  the  weight  selections  for  die  disturbance,  it  is  logical  to  utilize 
the  first  mediod  described  above.  If  die  system  is  assumed  to  be  that  of  a  typical 
aircraft,  a  disturbance  would  be  a  low  frequency  event  with  varying  magnitude.  With 
the  noise  signal  itself  assumed  to  be  white  noise,  the  disturbance  weight  will  be 
approximated  by  a  low  pass  filter  with  a  pole  at  s  =  -2.  For  comparison  purposes,  the 
magnitudes  of  this  weight  were  chosen  somewhat  arbitrarily,  corresponding  to  a  small, 
medium,  and  large  disturbance.  The  disturbance  weights  are  given  by 

„  0.0002  20  200  -X 

- -  W-'l 

s  +  2  s  +  2  s  +  2 

and  their  Bode  magnitude  plots  are  shown  in  Figure  S-3.  For  clarity,  however,  the 
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Figure  5-3.  Bode  Magnitude  Plot  of  Disturbance  Weights 


weights  will  be  described  on  example  plots  by  their  DC  gain,  which  are  -80  dB,  20 
dB,  and  40  dB,  respectively. 

The  measuremoit  noise  weights  utilize  the  same  logic  described  for  the 
disturbance  weights,  only  in  this  case,  measuremoit  noise  is  a  high  frequency 
phenomoion  with  varying  magnitude.  The  measurement  noise  weights  will  be 
q)proximated  by  a  high  pass  Alter  with  a  zero  at  -0.01  and  a  pole  at  -10.  The 
magnitudes  of  diis  weight  were  chosen  to  reflect  what  would  be  conadered  a  "small" 
measurmoit  noise  level  versus  a  "large"  noise  level.  The  measurement  noise  weights 
are  plotted  in  Figure  5-4  and  are  given  by 
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Figure  5-4.  Bode  Magnitude  Plot  of  Measurement  Noise  Weights 

^  ^  0.001(5+0.01)  20(5+0.01)  ^5  8^ 

On  example  plots,  these  two  noise  weights  will  be  referred  to  as  either  a  "small"  norse 
or  "large"  noise  to  reflect  the  magnitude  of  the  noise  weight 

The  only  remaining  dynamic  weight  for  this  ^stem  is  W^,  the  tracking  weight 
This  weight  is  of  the  type  that  should  be  chosen  r*:  a  penalty  on  the  signal  over  the 
frequency  range  of  interest.  As  the  name  implies,  the  tracking  weight  penalizes 
deviations  of  the  ou^ut  signal,  y,  from  the  reference  command,  r^.  Therefore,  in  order 
to  tqrply  this  penalty  over  the  relevant  frequency  range,  the  weight  should  be  shaped 
like  a  low  pass  filter  with  a  low  frequency  pole.  There  is  another  factor  to  consider. 
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however.  If  the  tracking  weight  is  set  to  be  a  low  pass  filter  (a  oHistant  gain  divided 
by  a  single  pole  near  the  origin),  thoi  the  weight  actually  resembles  an  integrator  over 
die  firequency  range  whoe  commands  come  in.  If  die  tracking  weight  resembles  an 
integrator,  then  the  actual  error  signal  will  reflect  the  inverse  of  the  integrator  and  die 
^stem  will  exhibit  type  I  behavior.  Likewise,  if  die  tracking  weight  is  set  to  be  a  low 
pass  filter  with  two  poles  near  the  origin,  then  the  ^ston  will  exhibit  type  n  behavior; 
and  so  on.  This  ^stem  type  bdiavior  will  be  examined  more  closely  in  an  example  in 
a  later  chiqiter.  Since  most  of  the  illustrative  examples  in  the  coming  cluqiters  concern 
themselves  only  with  step  re^onses,  die  tracking  weight  used  in  those  examples  is 


W, 


10 


s  *  0.01 


(5.9) 


Its  Bode  plot  is  shown  in  Figure  5-5. 


Figure  5-5.  Bode  Magnitude  Plot  of  Tracking  Weight 


As  mentioned  previously,  it  is  not  a  control  objective  to  minimize  the  2-norm 
of  the  to  y  transfer  function.  Thus,  it  would  not  make  sense  to  place  a  weight  on 
For  simplicity,  was  chosen  to  be  1.0.  Control  usage  is  often  a  limiting  factor 
in  aircraft  control  problems,  since  actuators  and  control  surfaces  have  limited  range 
and  speed.  Since  the  primary  purpose  of  diis  thesis  is  to  show  trends  in  the  speed  and 
robustness  of  1-DOF  and  2-DOF  systems  as  disturbances  and  measurement  noises  are 
varied,  the  control  usage  wei^t  was  set  to  a  value  which  allows  these  parameter 
variations  to  be  seen  easily  --  ITj  "  0  001. 

It  also  bears  mentioning  here  diat  the  weight  choices  described  previously  in 
this  section  are  applicable  to  the  multivariable  case.  The  primary  difference,  of 
course,  is  that  in  die  MIMO  case,  die  weights  represent  transfer  function  matrices 
instead  of  individual  transfer  functions.  While  it  is  true  that  it  is  often  desirable  to 
tailor  the  weights  to  each  individual  element  in  a  vector  signal,  it  is  not  necessary  for 
observing  trends.  For  the  purposes  of  this  diesis,  the  SISO  weights  will  be  multiplied 
by  an  appropriate  sized  identity  matrix  in  order  to  come  up  widi  a  MIMO  weight. 

5.4  Evaluation  Model 

Once  ail  of  the  weights  have  been  selected,  an  H2  algorithm  is  utilized  to 
calculate  the  optimal  compensator  K2g^  as  described  in  Chapter  4.  Since  the 
controller  is  designed  while  taking  all  of  the  weights  into  consideration,  it  would 
follow  that  if  a  simulation  were  set  up  that  incorporated  disturbances  and  measurement 
noises  of  intensity  and  frequency  similar  to  the  weights  used  in  the  design  of  the 
compensator,  then  the  system  with  the  newly  designed  controller  would  perform  quite 
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well.  For  the  purposes  of  this  Aesis,  however,  it  is  of  greater  interest  to  observe  the 
effects  of  weight  changes  on  both  1-DOF  and  2-DOF  systems  using  a  model  that 
purely  tests  the  performance  of  die  controller. 

The  evaluation  model  used  throughout  diis  thesis  is  simply  an  unweighted 
version  of  the  basic  feedback  setup  with  no  noises  or  "errors"  included.  The  block 
diagrams  of  the  evaluation  models  for  the  1-DOF  and  2-DOF  cases  are  shown  in 
Figures  S-6  and  S-7,  respectively.  These  may  not  be  "real  world"  simulation  models 
with  accurate  noises,  but  they  do  provide  a  consistent  basis  for  comparison  of  different 
controllers.  The  purpose  of  these  evaluation  models  is  not  to  evaluate  the  "goodness" 
of  an  H2  design;  rather,  it  is  to  examine  the  basic  differences  between  1-DOF  and  2- 
DOF  controllers. 

The  next  two  chapters  of  this  diesis  illustrate  the  concepts  described  thus  far 
through  examples  that  highlight  the  differences  between  the  1-DOF  and  2-DOF  design 
methodologies.  All  of  the  weights  given  in  this  chapter  will  be  used  in  each  problem 
to  maintain  consistency  for  comparison  purposes. 


Figure  5-6.  1-DOF  Evaluation  Model  Block  Diagram 
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VL  SISO  Examples 


To  best  illustrate  the  differences  between  1-DOF  and  2-DOF  control  systems,  it 
is  perii^s  best  to  explore  a  simple  example  that  is  representative  of  the  types  of 
problems  encountered  in  control  system  design.  In  this  chq>ter,  the  plant  consists  of  a 
transfer  function  widi  one  zero  and  two  poles.  This  example  can  be  manipulated  to 
produce  a  stable,  minimum  phase  system,  an  unstable,  minimum  phase  system,  a 
stable,  nonminimum  phase  ^stem,  and  an  unstable,  nonminimum  phase  system.  The 
transfer  functions  used  in  the  examples  are  as  follows: 


G(s) 

5+2 

s 

stable,  minimum  phase 

(6.1) 

(s  +  l)(s  +5) 

G(s)  ■■ 

5+2 

s 

unstable,  minimum  phase 

(6.2) 

(s  -  l)(s  +  5) 

G(s)  = 

s  -  2 

stable,  nonminimum  phase 

(6.3) 

(s  +  l)(s  +  5) 

G(s)  « 

s  -  2 

unstable,  nonminimum  phase 

(6.4) 

(s  -  l)(s  +  5) 

Although  these  transfer  functions  initially  appear  to  be  quite  simple,  it  turns  out  that  as 
far  as  general  trends  are  concerned  (regarding  the  effects  of  stable  vs.  unstable  poles 
and  minimum  phase  vs.  nonminimum  phase  zeros),  they  are  representative  of  the 
behavior  of  each  class  of  transfer  function. 

For  each  example  type,  there  are  many  gr^hs  of  loop  shapes  and  time 
responses  to  observe  in  order  to  make  comparisons  between  the  1-DOF  and  2<DOF 
model,  so  it  would  be  prudent  here  to  explain  the  order  in  vsdiich  the  gr^hs  will  be 
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presented.  In  each  class  of  problem,  die  1-DOF  plots  will  be  presented  first,  with 
plots  of  saisidvity,  complimentary  sensitivity,  loop  transfer  function,  and  step 
responses,  followed  by  the  same  plots  for  the  2-DOF  case.  Since  the  objectives  of  diis 
thesis  include  observing  similarities  and  differences  between  die  two  design 
methodologies  as  plant  disturbances  and  measurement  noises  are  varied,  each 
individual  plot  will  contain  three  curves,  corresponding  to  the  small,  medium,  and 
large  disturbances  defined  in  equation  (S.7).  As  mentioned  in  Chapter  V,  each  curve 
will  be  identified  by  the  appropriate  DC  gain  of  the  disturbance  (-80,  20,  and  40  dB, 
respectively).  Also,  the  graphs  will  appear  in  pairs,  with  "small"  measurement  noise 
plots  (as  defined  in  equation  (S.8))  being  on  the  top  of  the  page  and  "large" 
measurement  noise  plots  on  the  bottom. 

The  first  four  sections  of  this  chapter  will  concern  themselves  with  examining 
die  frequency  response  plots  and  step  responses  of  each  of  the  four  classes  of  example, 
making  relevant  comparisons  and  contrasts  i^ere  appropriate.  The  fifdi  section  of  the 
chapter  will  explore  the  effects  of  different  types  of  inputs  on  die  time  responses  of 
one  of  die  examples.  Specifically,  die  chosen  plant  will  be  given  step,  ramp,  and 
parabolic  inputs,  and  the  1-DOF  and  2-DOF  time  responses  will  be  compared.  The 
final  section  of  die  chtqiter  will  examine  the  possible  tracking  advantages  of  using  a  2- 
DOF  compensator  \^en  the  system  is  subjected  to  measurement  noise  which  contains 
low  frequency  content. 

6.1  Stable,  Minimum  Phase  System 

The  stable,  minimum  phase  example  serves  as  a  baseline  for  comparison  to  all 
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of  tiie  otfier  types  of  examples,  but  it  is  still  instructive  to  analyze  all  of  the  plots  in 
order  to  reiterate  the  concepts  from  Chapters  n  and  m. 

6.1.1  l-DOF  Case.  For  the  basic  1-DOF  case,  all  of  the  loop  shape  plots 
should  closely  resemble  die  plots  of  "good"  loop  shapes  in  Chapter  n.  Beginning  with 
die  sensitivity  function.  Figure  6-1  shows  that  is  small  at  low  frequracy  and  levels 
out  at  0  dB  at  higher  frequencies.  As  die  disturbance  becomes  larger,  \GK\  also 
becomes  larger  at  low  frequency.  In  the  large  measurement  noise  case  in  Figure  6-2, 
the  shape  of  |•S|  is  exacdy  the  same,  but  it  does  not  have  as  small  a  magnitude  at  low 
frequency.  This  is  expected  due  to  the  fact  diat  with  a  large  measurement  noise,  \GK\ 
must  be  smaller  at  high  frequency  to  reject  it,  and  the  overall  magnitude  of  GK  is 
lower  over  all  frequencies. 

The  complimentary  sensitivity  plots  reflect  that  when  \GK\  is  large  at  low 
frequency,  17]  iqiproaches  unity,  or  0  dB,  and  rolls  off  at  high  frequency.  As  the 
disturbance  becomes  larger,  |G^  becomes  larger  at  low  frequency,  and  the  bandwidth 
of  the  system  increases.  A  comparison  of  Figures  6-3  and  6-4  verifies  the  effect  of  a 
large  measurement  noise  on  the  system  —  |7]  rolls  off  sooner  when  the  noise  is  large 
because  \GK\  is  smaller  at  high  frequency. 

Figures  6-5  and  6-6  illustrate  the  shape  of  a  "good"  loop  transfer  function,  with 
high  gain  at  low  frequency  and  low  gain  at  high  frequency.  Once  again,  \GK\ 
increases  as  the  disturbance  gets  larger.  The  effect  of  die  large  measurement  noise  is 
to  shift  die  magnitude  of  \GK\  downward  for  each  disturbance. 

Figures  6-7  and  6-8  feature  the  step  responses  of  the  1-DOF  system,  utilizing 
the  evaluation  model  from  Chapter  V,  with  no  plant  disturbances  or  measurement 
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noises  acting  on  the  system.  As  the  magnitude  of  the  assumed  disturbance  increases, 
so  does  die  magnitude  of  GK  in  the  resulting  controller  at  low  frequwcy,  which  makes 
the  steady  state  error  smaller.  Differences  in  steady  state  error  are  not  visible  in  the 
figures  because  they  are  so  small,  however.  The  bandwiddi  frequency  also  increases 
as  \GK\  gets  larger,  vdiich  has  a  much  more  noticeable  effect  on  the  response  by 
making  the  rise  times  much  faster.  In  the  large  measurement  noise  case,  the  same 
trends  are  present,  but  since  \GK\  is  not  as  large  at  low  frequency,  and  the  bandwiddi 
frequency  is  lower  for  each  disturbance,  die  setding  times  are  much  longer.  When  the 
disturbance  is  large  enough,  the  system  steadies  out  in  less  than  a  second,  but  the  large 
noise  manifests  itself  by  making  the  initial  overshoot  larger  than  that  of  the 
corresponding  small  noise  case. 
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Magniinde  (dB)  ^  Magulnde  (dB) 


Igure  6-1.  1-DOF  Bode  Magnitude  Plot  of  151  for  Stable,  Minimum  Phase 
System  Example  (Small  Noise  Case) 


Figure  6-2.  1-DOF  Bode  Magnitude  Plot  of  151  for  Stable,  Minimum  Riase 
System  Example  (Large  Noise  Case) 
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Frequotcy  (nd/m) 


Figure  6*3.  l-DOF  Bode  Magnitude  Plot  of  171  for  Stable,  Minimum 
Phase  System  Example  (Small  Noise  Case) 


Frequency  (lud/sec) 

Figure  6-4.  1-DOF  Bode  Magnitude  Plot  of  171  for  Stable,  Minimum 
Phase  System  Example  (Large  Noise  Case) 
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Magnitude  (dB)  Magnitude  (dB) 


Figure  <•?.  I'DOF  Step  Re^Kmses  for  Stable,  Minimum  Phase  Systran 
Example  (Small  Noise  Case) 


Figure  6>8.  l-DOF  Step  Responses  for  Stable,  Minimum  I%ase  System 
Example  (Large  Noise  Case) 
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6,1.2  2~DOF  Case.  The  2-DOF  case  also  serves  as  a  baseline  for  the  otfier 
types  of  2>DOF  examples  in  that  the  loop  shape  plots  resemble  those  of  "good"  loop 
shapes  as  described  in  Chapter  m.  Figures  6-9  and  6-10  show  magnitude  plots  of  the 
sensitivity  function.  In  both  the  small  and  large  noise  cases,  |5|  exhibits  the 
characteristic  "good"  shape  of  being  small  at  low  frequency  and  leveling  out  at  0  dB  at 
high  frequency.  As  the  assumed  plant  disturbance  level  is  increased,  the  magnitude  of 
GK  is  increased  at  low  frequency;  thus,  the  magnitude  of  |5)  decreases  at  low 
frequency.  In  the  large  noise-small  disturbance  case,  |<S]  is  at  0  dB  at  virtually  all 
frequencies  because  the  large  noise  f  tees  dovsm  the  size  of  \GK^\  over  all  frequencies, 
and  the  disturbance  is  not  large  enough  to  offset  this. 

Figures  6-11  and  6-12  show  plots  of  die  2-DOF  complimentary  sensitivity 
function.  These  plots  also  exhibit  the  characteristics  of  a  good  loop  shape,  in  that 
most  of  the  plots  are  at  0  dB  at  low  freque;.  .:y  and  roll  off  at  some  bandwidth 
frequency.  In  the  small  noise  case,  the  disturbance  is  always  large  enough  for  jT]  to  be 
0  dB  at  low  frequency.  As  the  disturbance  gets  larger,  the  bandwidth  frequency 
increases.  In  the  large  noise  case,  however,  the  size  of  the  noise  holds  down  \GK^\ 
since  the  disturbance  is  small,  and  |21  is  not  even  guaranteed  to  be  at  0  dB  at  low 
frequency.  In  fact,  when  the  disturbance  is  at  its  smallest,  the  DC  gain  of  |7]  is  nearly 
-100  dB.  The  roll-off  frequencies  for  all  of  the  disturbances  are  significantly  lower  for 
the  large  noise  case  than  those  for  the  small  noise  case. 

Figures  6-13  and  6-14  show  plots  of  \GK^\  and  \GK^  for  the  small  noise  case. 
The  plots  of  1GK,1  reflect  the  trends  observed  in  the  plots  of  151  and  \T\.  with  high  gain 
at  low  frequency  and  low  gain  at  high  frequency.  As  the  size  of  the  disturbance 
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increases,  Ae  overall  size  of  IG/^il  also  increases.  For  Ae  two  larger  Asturbances, 
\GKJt  viiiiially  matches  \GKil  but  in  Ae  small  Asturbance  case, 
compensate  for  thft  effects  of  Ae  measurement  noise  by  flattening  out  at  0  dB  until 
around  10  rad/sec  before  rolling  off.  In  Ae  large  noise  case  (Figures  6-15  and  6-16), 
Ae  overall  magnitudes  of  IG/Tjl  decrease  much  more  rapiAy  wiA  decreasing 
Asturbances  Aan  in  Ae  1-DOF  case,  especially  in  Ae  smallest  Asturbance  case,  where 
IG/Tjl  has  a  DC  gain  of  -100  dB.  Since  Ae  overall  system  bandwidA  of  Ae  2-DOF 
system  is  about  10  rad/sec,  (GATjl  shapes  itself  to  roll  off  at  Ais  frequency,  even  if 
IG/Tjl  falls  below  0  dB  before  Ais  frequency.  If  Ae  Asturbance  is  large  enough 
(relative  to  Ae  noise),  it  means  Aat  IG/TJ  will  probably  match  Ae  shape  of  \GK^[,  if 
Ae  Asturbance  :  sMall  enough,  Aen  IG/T^i  probably  have  to  stay  0  dB  until  it 
reaches  10  rad/sec.  This  trend  is  repeated  here  because  it  will  recur  in  all  of  Ae 
remaining  examples  in  Ais  Aesis. 

Figure  6-17  shows  a  plot  of  Ae  system  transfer  function  from  to  y,  given  by 
which  (from  Chapter  ID)  is  Ae  2-DOF  eqmvalent  of 'T)  for  Ae  1-DOF  case 
insofar  as  evaluating  Ae  frequency  range  over  vAich  Ae  system  will  track  commands. 
It  shows  Aat  Ae  2-DOF  system  bandwidA  for  command  following  is  indeed  about  10 
raAsec. 

The  final  graph  for  Ae  stable,  minimum  phase  case.  Figure  6- 1 8,  is  Ae  step 
response  for  Ae  2-DOF  system.  It  is  immeAately  apparent  Aat  Aere  is  only  one  line 
on  Ae  graph,  as  Aere  always  is  for  2-DOF  controllers.  WiA  Ae  flexibility  of  two 
compensators,  Ae  2-DOF  system  is  able  to  produce  Ae  optimal  time  response, 
regardless  of  Ae  size  of  Ae  Asturbance  or  measurement  noise.  In  Ais  case,  Ae 
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response  exhibits  Ae  shq)e  of  a  critically  damped  system,  with  the  settling  time  being 
as  fast  as  possible  without  any  overshoot . 

6.1.3  Comparison  Between  1-DO F  and  2-DO F  Cases.  The  differences 
between  the  1-DOF  and  2-DOF  cases  for  this  first  basic  example  are  centered  around 
the  different  requirements  for  good  loop  shapes,  as  described  in  Chapters  II  and  in. 
The  2-DOF  sensitivity  function  shows  the  same  trends  as  in  the  1-DOF  example  in  the 
small  noise  case,  except  that  for  the  smallest  disturbance,  it  is  not  as  small  at  low 
frequency  as  it  is  in  the  1-DOF  case.  The  reason  for  |5|  not  being  as  small  at  low 
frequency  in  the  2-DOF  case  stems  from  the  idea  that  \GK^\  is  only  as  large  as  it  needs 
to  be  to  reject  the  l^}propriate  sized  disturbance  at  low  frequency.  In  both  the  small 
and  large  noise  cases,  |5]  is  at  least  as  small  at  low  frequency  as  the  corresponding  1- 
DOF  case  vdien  the  disturbance  is  at  its  largest,  reflecting  the  fact  that  \GK^\  must  be 
sufficiently  large  at  low  frequency  to  reject  a  large  disturbance.  The  2-DOF  Bode 
magnitude  plot  of  |7|  differs  from  that  of  the  1-DOF  case  in  that  it  is  not  even 
guaranteed  to  be  at  0  dB  at  low  frequency  vriien  the  measurement  noise  level  is  large, 
due  to  the  fact  that  IG/Til  is  held  down  over  all  frequencies.  Since  the  2-DOF  loop 
transfer  function,  GK^,  does  not  have  the  same  restrictions  on  command  following  that 
the  1-DOF  loop  transfer  function  does,  the  plots  of  |GA^||  show  the  variability  in  the 
magnitudes,  depending  on  the  size  of  the  assumed  disturbance. 

The  rise  time  and  settling  time  of  the  2-DOF  step  response  are  not  as  fast  as 
some  of  the  1-DOF  responses  in  the  small  noise  case,  but  they  are  faster  than  those  of 
the  large  noise  case  (note  the  time  scale  differences  in  the  figures).  The  advantage  of 
the  2-DOF  controller  is  that  although  the  response  is  not  as  fast  as  some  1-DOF 
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special  cases,  it  always  produces  the  same  response  (for  a  given  tracking  weight),  no 
matter  >^at  the  size  of  the  disturbance  or  measurement  noise  is.  On  the  other  hand,  a 
large  measurement  noise  has  a  detrimental  effect  on  the  1-DOF  responses.  It  should 
be  noted  diat  the  shq>e  of  the  2-DOF  step  response  can  be  altered  by  changing  the 
tracking  weight.  For  instance,  if  the  gain  of  the  tracking  weight  is  increased  so  dtat  its 
magnitude  plot  crosses  the  0  dB  line  at  a  higher  frequency,  the  resulting  step  response 
will  have  a  faster  rise  time,  but  more  overshoot.  Most  of  the  examples  in  this  thesis, 
however,  have  been  completed  using  die  tracking  weight  mentioned  in  Chtqiter  V, 
since  this  weight  produces  step  responses  in  which  the  differences  between  the  1-DOF 
and  2-DOF  case  can  readily  be  seen. 
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Figure  6>9.  2-DOF  Bode  Magnitude  Plot  of  ISI  for  Stable,  Minimum 
niase  Syston  Example  (Small  Noise  Case) 


Figure  6>10.  2-DOF  Bode  Magnitude  Plot  of  LSI  for  Stable,  Minimum 
Phase  System  Example  (Large  Noise  Case) 
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Figure  6-11 ,  2-DOF  Bode  Magnitude  Plot  of  171  for  Stable,  Minimum 
Phase  Syston  Example  (Small  Noise  Case) 


Figure  6-12.  2-DOF  Bode  Magnihide  Plot  of  171  for  Stable,  Minimum 
Phase  System  Example  (Large  Noise  Case) 
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Figure  6-13.  2-DOF  Bode  Magnitude  Plot  of  \GK^  for  Stable,  Minimum 
Phase  System  Example  (Small  Noise  Case) 


Hrequeocy  (rad/aec) 


Figure  6-14.  2-DOF  Bode  Magnitude  Plot  of  foi*  Stable,  Minimum 
Phase  System  Example  (Small  Noise  Case) 
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Figure  6«15.  2-DOF  Bode  Magnitude  Plot  of  for  Stable,  Minimum 
Hiase  Syston  Example  (Large  Noise  Case) 


Figure  6-16.  2-DOF  Bode  Magnitude  Plot  of  IGATjl  for  Stable,  Minimum 
niase  System  Example  (Large  Noise  (^ase) 
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Frequoicy  (radl/sec) 

Figure  6-17.  2-DOF  Bode  Magnitude  Plot  of  ISGiiTjl  for  Stable,  Minimum 
niase  Syston  Example  (any  size  noise) 


Figure  6-18.  2-DOF  Step  Response  for  Stable,  Minimum  I^ase  System 
Example  (any  size  noise) 
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6.2  Unstable,  Minimum  Phase  System 

62.1  I’DOF  Case.  In  this  case,  the  only  difference  in  the  plant  transfer 
function  is  that  now  there  is  an  unstable  pole  at  s  =  +1.  Figures  6-19  and  6-20  show 
plots  of  the  1-DOF  sensitivity  function,  |iS|.  In  both  the  small  and  large  noise  cases, 
there  is  no  significant  change  from  the  stable,  minin  .un  phase  example  in  trends  or  in 
the  actual  magnitudes  of  the  function.  Figures  6-21  and  6-22  show  plots  of  the 
complimentaiy  sensitivity  function,  |7].  Again,  the  same  trends  are  present  as  in  die 
stable,  minimum  phase  example,  except  that  the  function  exhibits  a  small  increase  in 
magnitude  before  rolling  off,  especially  in  the  large  noise  case.  This  "bubble" 
indicates  deteriorating  stability  margins,  which  should  be  expected  since  the  plant  is 
unstable  to  begin  with,  and  the  stability  margins  will  not  be  as  good  as  they  would  be 
with  a  stable  plant.  The  plots  of  \GK\  are  shown  in  Figures  6-23  and  6-24.  They  are 
virtually  identical  to  the  plots  of  \GK[  in  the  stable,  minimum  phase  case. 

The  real  differences  between  the  stable,  minimum  phase  and  the  unstable, 
minimum  phase  systems  are  evident  in  the  step  responses,  shown  in  Figures  6-2S  and 
6-26.  In  the  small  noise-large  disturbance  case,  the  responses  are  nearly  identical,  but 
when  the  disturbance  gets  smaller,  the  responses  begin  to  settle  out  in  the 
neighborhood  of  two  seconds  in  the  unstable  system  instead  of  one  second  in  the 
stable  system.  In  the  large  noise  case,  the  setding  times  are  generally  much  slower, 
ranging  from  just  under  a  second  for  the  large  dishirbance  case  to  over  10  seconds  for 
die  small  disturbance  case.  As  the  assumed  plant  disturbance  gets  smaller,  the  settling 
time  produced  by  the  resulting  compensator  gets  longer,  and  the  peak  overshoot  gets 
larger. 


6-18 


Mi^itnde  (dB)  Mi^goiliide  (dB) 


r 


Figure  6*19.  1-DOF  Bode  Magnitude  Plot  of  LSI  for  Unstable,  Minimum 
Phase  System  Example  (Small  Noise  C^ase) 


Figure  6>20.  1-DOF  Bode  Magnitude  Plot  of  LSI  for  Unstable.  Minimum 
Phase  System  Example  (Large  Noise  Case) 
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Figure  6-21.  1-DOF  Bode  Magnitude  Plot  of  171  for  Unstable,  Minimiun 
Phase  Syston  Example  (Small  Noise  Case) 


Figure  6«22.  1-DOF  Bode  Magnitude  Plot  of  171  for  Unstable,  Minimum 
Phase  System  Example  (Large  Noise  Case) 
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Figure  6*25.  l*DOF  Step  Reqionses  for  Unstable,  Minimum  Phase 
System  Example  (Small  Noise  Case) 


Figure  6-26.  1-DOF  Step  Responses  for  Unstable,  Minimum  Phase 
System  Example  (Large  Noise  Case) 


6.2.2  2~DOF  Case.  Figures  6-27  and  6-28  show  plots  of  the  2-DOF 
sensitivity  function.  They  are  virtually  identical  to  the  equivalent  plots  in  the  stable, 
minimum  phase  system.  Figures  6-29  and  6-30  show  plots  of  the  complimentaiy 
sensitivity  function.  These  plots  show  a  slight  increase  in  magnitude  just  before 
rolling  off  in  the  small  noise  case,  indicating  the  slight  deterioration  of  the  stability 
margins  compared  to  the  stable  case.  In  the  large  noise-small  disturbance  case,  the 
DC  gain  of  |71  is  slightly  above  0  dB  instead  of  there  being  only  a  slight  hump,  again 
indicating  worse  margins.  In  Figures  6-31  and  6-32,  the  plots  of  \GKi\  and  \GK^  for 
the  small  noise  case  are  largely  unchanged  from  the  stable  example,  as  are  the  plots 
for  the  large  noise  case,  in  Figures  6-33  and  6-34.  Once  again,  \GK.^  shapes  itself  to 
ensure  that  it  rolls  off  below  0  dB  at  least  10  rad/sec.  The  plot  of  \SGK^  is  identical 
to  Figure  6-17  and  shows  the  2-DOF  system  bandwidth  for  tracking  command  inputs. 
The  2-DOF  step  response  is  also  identical  to  that  of  the  2-DOF  stable  example  in 
Figure  6-18. 

The  primary  difference  between  the  stable,  minimum  phase  system  and  die 
unstable,  minimum  phase  case  is  the  shape  of  the  complimentary  sensitivity  plot.  In 
the  unstable  case,  the  function  17]  shows  an  increase  in  magnitude  before  rolling  off, 
indicating  a  worsening  of  the  stability  margins,  which  is  to  be  expected  vdien  going 
from  a  stable  to  an  unstable  plant. 

6.2.3  Comparison  Between  1-DOF  and  2~DOF  Cases.  All  of  the  loop  shape 
plots  for  the  1-DOF  and  2-DOF  systems  are  virtually  the  same  as  they  were  in  the 
stable  system  example,  except  for  the  complimentary  sensitivity  function  plots.  The  2- 
DOF  sensitivity  function  plots  again  show  that  since  IG/Tjl  is  not  required  to  be  large 
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at  low  frequency,  15]  is  larger  at  low  frequency  than  the  corresponding  1-DOF  case 
^en  the  disturbance  is  small  relative  to  the  measurement  noise.  In  both  the  1>D0F 
and  2-DOF  cases,  |7]  exhibits  a  rise  in  magnitude  at  some  point  before  rolling  off, 
indicating  a  decrease  in  stability  margins.  In  the  small  disturbance-large  noise  case, 
the  DC  gain  of  |7]  is  above  0  dB  in  the  2-DOF  case,  vdiich  never  occurs  in  the  1-DOF 
case,  due  to  the  requirement  for  the  1-DOF  loop  gain  to  be  large  at  low  frequency  (for 
both  disturbance  rejection  and  tracking). 

Overall,  an  unstable  pole  has  virtually  no  effect  on  the  tracking  ability  of  die  2- 
DOF  system,  while  it  has  noticeable  effects  on  the  performance  of  the  1-DOF  system. 
The  1-DOF  step  responses  show  that  an  unstable  pole  generally  increases  the  amount 
of  overshoot.  In  the  2-DOF  case,  |G^2l  shtqied  in  such  a  way  that  perfect  tracking 
is  maintained  (from  the  plot  of  |5Gif2l)  over  a  larger  frequency  range  than  the  1-DOF 
case. 
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Figure  6*27.  2-DOF  Bode  Magnitude  Plot  of  LSI  for  Unstable,  Minimum 
Phase  Syston  Example  (Small  Noise  Case) 
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Figure  6-28.  2-DOF  Bode  Magnitude  Plot  of  LSI  for  Unstable,  Minimum 
I%ase  System  Example  (Large  Noise  Case) 
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Figure  6*29.  2-DOF  Bode  Magnitude  Plot  of  171  for  Unstable,  Minimum 
Phase  System  Example  (Small  Noise  Case) 


Figure  6-30.  2-DOF  Bode  Magnitude  Plot  of  171  for  Unstable,  Minimum 
I%ase  System  Example  (Large  Noise  Case) 
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igure  6>31. 


2-DOF  Bode  Magnitude  Plot  of  IGiCil  for  Unstable,  Minimum 
I%ase  System  Example  (Small  Noise  (^ase) 


Figure  6>32.  2-DOF  Bode  Magnitude  Plot  of  KJKJi  for  Unstable,  Minimum 
I%ase  System  Example  (Small  Noise  Case) 
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igure  6-33.  2-DOF  Bode  Magnitude  Plot  of  IGAT,!  for  Unstable,  Minimum 
Phase  System  Example  (Large  Noise  Case) 
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Figure  6-34.  2-DOF  Bode  Magnitude  Plot  of  IGATjl  Unstable,  Minimum 
Phase  System  Example  (Large  Noise  Case) 
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6.3  Stable,  Nonminimum  Phase  System 

The  next  two  examples  presented  repres^t  nonminimum  phase  systems,  vdiich 
often  pose  challenging  control  problems.  Using  feedback  control,  poles  in  undesirable 
locations  can  be  moved  to  more  favorable  closed  loop  locations,  but  undesirable  zeros 
are  a  burden  which  cannot  be  directly  remedied  —  zeros  cannot  be  re-assigned  to 
different  locations. 

6.3.1  1-DOF  Case.  Figures  6-3S  and  6-36  show  plots  of  the  sensitivity 
function,  jiS).  The  most  noticeable  characteristic  of  these  plots  is  the  large  amount  of 
overshoot  exhibited  by  the  functions.  The  existence  of  a  nonminimum  phase  zero  has 
the  effect  of  attracting  stable  poles  toward  the  right-half  plane  (on  a  root  locus  plot)  as 
gain  is  increased;  therefore,  the  stability  margins  are  expected  to  be  lower  for  this 
example  dian  either  of  the  two  minimum  phase  examples  described  thus  far.  The  plots 
of  |7]  in  Figures  6-37  and  6-38  reiterate  the  decreased  stability  margins  in  that  they 
also  show  a  significant  rise  in  magnitude  before  rolling  off  at  the  appropriate 
bandwidth  frequency.  The  amount  of  overshoot  in  |5|  and  |71  is  a  bit  more  pronounced 
in  the  small  noise  case  dian  it  is  in  the  large  noise  case.  This  is  because  in  the  small 
noise  case,  die  size  of  \GK\  is  permitted  to  be  larger  overall,  thereby  placing  the  closed 
loop  poles  in  a  position  where  the  response  is  faster,  but  some  pole(s)  may  be  closer 
to  instability. 

Figures  6-39  and  6-40  show  plots  of  die  loop  shape,  \GK\.  Both  plots  still 
exhibit  high  gain  at  low  frequency,  but  the  effect  of  the  nonminimum  phase  zero  is 
evident  in  their  sh{q)es.  Since  zeros  cannot  be  moved  under  feedback  control,  the  next 
best  manner  in  vdiich  to  handle  the  situation  is  to  minimize  the  amount  of  energy 
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placed  at  die  frequency  of  die  nonminimum  phase  zero.  This  is  shown  in  the  two 
figures  by  die  fact  that  none  of  the  lines  is  idiove  0  dB  past  2  rad/sec,  which 
corresponds  to  the  zero  at  s  =  +2.  In  fact,  that  particular  frequency  effectively 
presents  a  "barrier"  under  which  the  magnitude  of  GK  must  pass  in  order  to  produce 
stable  closed  loop  systems.  In  the  small  noise  case,  \GK\  normally  has  higher  gain  at 
low  frequency  and  does  not  cross  the  0  dB  line  until  at  least  10  rad/sec.  For  bodi 
noise  levels  with  a  nom.iiuimum  phase  zero,  however,  \GK\  passes  under  the  2  rad/sec 
barrier  and  then  stays  just  below  0  dB  out  to  frequencies  beyond  100  rad/sec  before 
rolling  off,  thus  trying  to  maintain  the  highest  amount  of  energy  possible  over  all 
frequencies. 

Figures  6-41  and  6-42  show  the  step  responses  for  the  stable,  nonminimum 
phase  system.  As  expected,  the  initial  tendency  of  nonminimum  phase  systems  is  to 
initially  "go  the  wrong  way."  Once  again,  the  general  trends  are  as  before  -  the 
setding  times  get  faster  as  the  disturbance  gets  bigger  in  both  the  large  and  small  noise 
cases,  but  the  large  noise  case  produces  slower  setding  times  than  the  small  noise  case 
does.  In  the  large  disturbance  case,  the  settling  times  are  about  the  same  for  both 
noises,  but  the  initial  undershoot  is  much  larger  in  the  small  noise  case.  This  is  in 
keeping  with  the  fact  that  as  |G^  is  made  larger  in  a  nonminimum  phase  system,  the 
response  gets  faster,  but  some  pole(s)  may  move  closer  to  the  imaginary  axis,  causing 
more  oscillatory  motion.  When  the  disturbances  are  smaller  in  the  large  noise  case, 
there  is  not  as  much  undershoot,  but  the  settling  times  become  quite  a  bit  slower. 
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i'igurc  6-35.  1-DOF  Bode  Magnitude  Plot  of  151  for  Stable,  Nonminimum 
niase  System  Example  (Small  Noise  Case) 


Figure  6-36.  1-DOF  Bode  Magnitude  Plot  of  151  for  Stable,  Nonminimum 
Phase  System  Example  (Large  Noise  Case) 
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Figure  6>37.  l«DOF  Bode  Magnitude  Plot  of  171  for  Stable,  Nonminimum 
Phase  System  Example  (Small  Noise  Case) 


Figure  6-38.  1-DOF  Bode  Magnitude  Plot  of  171  for  Stable,  Nonminimum 
Phase  System  Example  (Large  Noise  Case) 
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igure  6-39. 


1-DOF  Bode  Magnitude  Plot  of  IGiin  for  Stable,  Nonminimum 
Phase  System  Example  (Small  Noise  Case) 


Figure  6-40.  1-DOF  Bode  Magnitude  Plot  of  iGi^l  for  Stable,  Nonminimum 
Phase  Syst^  Example  (Large  Noise  Case) 
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Figure  6>41.  1-DOF  Step  Refuses  for  Stable,  Nonminimum  I%ase 
System  Example  (Small  Noise  Case) 


Figure  6-42.  1-DOF  Step  Responses  for  Stable,  Nonminimum  Phase 
System  Example  (Large  Noise  Case) 
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6.3.2  2~DOF  Case.  The  2-DOF  sensitivity  functions  are  plotted  in  Figures  6- 
43  and  6-44.  The  general  shape  of  the  plots  is  the  same  as  that  of  "good"  loop  shapes 
(from  Chapter  m),  but  there  is  a  significant  amount  of  overshoot  in  both  the  small  and 
large  noise  cases.  The  functions  in  the  small  noise  case  have  a  bit  more  overshoot 
than  those  in  the  large  noise  case,  due  the  fact  that  the  loop  gain  is  higher  overall.  In 
the  large  noise-small  disturbance  case,  the  size  of  \GK^\  is  so  small  that  |5]  is  at  0  dB 
at  virtually  all  frequencies.  Figures  6-45  and  6-46  show  plots  of  the  complimentary 
sensitivity  function.  These  plots  also  indicate  the  decrease  in  stability  margins  by  the 
large  increase  in  magnitude  before  rolling  off.  The  other  trends  for  these  plots  are  in 
agreement  with  the  previous  examples. 

Figures  6-47  and  6-48  show  plots  of  \GK^\  and  jGATJ  for  the  small  noise  case. 

In  the  small  disturbance  case,  \GK^  exhibits  high  gain  at  low  frequency  and  low  gain 
at  high  frequency,  but  since  the  assumed  disturbance  level  is  so  low,  the  plot  is  not 
affected  by  the  barrier  at  2  rad/sec.  In  the  two  larger  disturbance  cases,  however, 

\GK^\  is  required  to  be  large  at  low  frequen^r,  but  it  must  also  pass  under  the  barrier 
at  2  rad/sec.  Figures  6-49  and  6-50  show  plots  of  |GAr,|  and  \GK2\  for  the  large  noise 
case.  The  effect  of  the  large  noise,  once  again,  is  to  drive  down  the  size  of  IG/Til, 
especially  in  the  small  distiutance  case  The  barrier  at  2  rad/sec  is  evident  in  the  plots 
of  both  IGATjl  and  IG/TJ-  Figure  6-51  shows  a  plot  of  \SGK'^,  indicating  the  frequency 
range  over  wdiicb  the  system  will  track  commands.  Figure  6-52  sho^s  the  2-DOF  step 
response. 

6.3.3  Comparison  Between  1-DOF  and  2-DOF  Cases.  In  keeping  with 
previous  trends,  the  plots  of  the  sensitivity  function  show  close  agreement  between  the 
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1-DOF  and  2-DOF  cases  when  the  measurement  noise  level  is  small.  In  the  large 
noise  case,  the  DC  gain  of  liS]  increases  more  r{^>i(lly  in  the  2-DOF  case  than  it  does  in 
the  1-DOF  case  as  the  size  of  the  assumed  disturbance  decreases.  The  plots  of  die 
complimentary  sensitivity  fimction  also  follow  previous  trends,  with  the  bandwidth 
frequency  increasing  as  the  assumed  disturbance  gets  larger  in  the  small  noise  case.  In 
the  large  noise  case,  the  same  trends  are  present,  but  in  the  2-DOF  large  noise-small 
disturbance  case,  die  DC  gain  of \J]  is  -100  dB.  In  the  small  noise  case,  the  plots  of 
\GKi\  are  very  similar  to  the  plots  of  |G/C|  in  the  1-DOF  case,  but  the  2-DOF  plot  for 
the  largest  disturbance  case  rolls  off  at  a  higher  frequency  than  the  corresponding  1- 
DOF  plot.  This  is  because  |G^i|  is  not  constrained  to  any  tracking  requirement,  and 
the  measurement  noise  is  not  large  enough  to  warrant  holding  down  its  size.  Since  the 
measurement  noise  is  so  small,  IG/Tjl  tends  to  dominate  the  system's  behavior;  thus, 
the  shqie  of  |G#r2l  is  more  or  less  the  same  as  that  of  |G^||.  It  should  be  noted  that 
both  |G^i|  and  |G^2l  ™ust  fall  beneath  the  barrier  at  2  rad/sec  to  accommodate  the 
zero  at  s  =  +2.  In  the  large  noise  case,  the  shape  of  \GKi\  and  \GK2\  follow  the  same 
pattern  that  they  do  in  the  other  classes  of  examples. 

The  2-DOF  step  response's  settling  time  of  about  three  seconds  is  not  as  fast  as 
some  of  the  1-DOF  responses  (in  the  small  noise  case),  but  it  does  not  have  as  much 
initial  undershoot  or  overshoot  of  the  final  value  as  they  do.  Once  again,  the  2-DOF 
system  produces  the  same  step  response,  regardless  of  the  level  of  the  plant 
disturbance  or  measurement  noise  level. 
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Figure  <>«43.  2<DOF  Bode  Magnitude  Plot  of  LSI  for  Stable,  Nonmmimuni 
I%ase  Systm  Example  (Small  Noise  Case) 


Frequency  (rad/sec) 

Figure  6>44.  2-DOF  Bode  Magnitude  Plot  of  151  for  Stable,  Nonmii  num 
Phase  System  Example  (Large  Noise  Case) 
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Figure  6*45.  2'DOF  Bode  Magnitude  Plot  of  171  for  Stable,  Nonminimum 
I%ase  System  Example  (Small  Noise  C^) 


Figure  6>46.  2-DOF  Bode  Magnitude  Plot  of  171  for  Stable,  Nonminimum 
I%ase  System  Example  (Large  Noise  Case) 
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Figure  6<47.  2-IX)F  Bode  Magnitude  Plot  of  IGiiT,!  for  Stable,  Nonminimom 
I%ase  System  Example  (Small  Noise  Case) 


Figure  6-48.  2-DOF  Bode  Magnitude  Plot  of  \GK^  for  Stable,  Nonminimom 
Phase  System  Example  (Small  Noise  Case) 


6-39 


Figure  6-49.  2-DOF  Bode  Magnitude  Plot  of  KjAT,!  for  Stable,  Nonminimum 
Phase  System  Example  (Large  Noise  Case) 


Figure  6-50.  2-DOF  Bode  Magnitude  Plot  of  IGiiTjl  for  Stable,  Nonminimum 
Phase  System  Example  (Large  Noise  Case) 
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Figure  6*51.  2-DOF  Bode  Magnibide  Plot  of  \SGKji  for  Stable, 

Nonminimum  I^iase  System  Example  (any  size  noise) 


Time  (sec) 

Figure  6-52.  2-DOF  Step  Re^onse  for  Stable,  Nonminimum  Phase 
System  Example  (any  size  noise) 
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6,4  Unstable f  Nonminimum  Phase  System 

6.4,1  1-DOF  Case.  Figures  6-53  and  6-54  show  plots  of  die  sensitivity 
function.  They  are  very  similar  to  the  stable,  nonminimum  phase  case,  but  the 
functions  do  not  get  as  small  at  low  frequency,  and  die  overshoot  is  much  more 
pronounced.  The  presence  of  right-half  plane  poles  and  zeros  erodes  the  stability 
margins  more  dian  in  any  odier  case,  and  the  size  of  |G^  is  held  down  overall  as  a 
result.  The  grtqihs  of  the  complimentaiy  sensitivity  function,  in  figures  6-55  and  6-56, 
add  further  evidence  of  low  stability  margins,  with  bigger  increases  in  magnitude 
before  rolling  off.  The  shiqie  of  \GK\  (Figures  6-57  and  6-58)  reiterates  the 
requirement  that  the  function  stay  below  0  dB  at  the  2  rad/sec  nonminimum  phase  zero 
frequency.  In  die  small  noise-large  disturbance  case,  the  function  "lips  up"  above  0 
dB  at  2000  rad/sec  in  an  effort  to  regain  some  overall  energy  before  rolling  off,  even 
though  it  still  avoids  violating  the  barrier  at  2  rad/sec. 

The  1-DOF  step  responses  are  shown  in  Figures  6-59  and  6-60.  These 
responses  exhibit  characteristics  of  both  unstable  and  nonminimum  phase  behavior,  as 
expected,  and  show  the  same  general  trends  as  the  stable,  nonminimum  phase  case.  In 
the  small  noise  case,  the  initial  undershoot  increases  as  the  level  of  die  assumed  plant 
disturbance  increases.  In  die  large  noise  case,  the  initial  undershoot  also  increases,  but 
the  subsequent  overshoot  of  the  final  value  decreases  slighdy,  as  the  assumed 
disturbance  level  increases. 
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igure  6-53.  1-DOF  Bode  Magnitude  Plot  of  151  for  Unstable,  Nonnunimum 
Phase  System  Example  (Small  Noise  Case) 


Figure  6*54.  1-DOF  Bode  Magnitude  Plot  of  151  for  Unstable,  Nonminimum 
Phase  System  Example  (Large  Noise  Case) 
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igure  6>S5. 


l.DOF  Bode  Magnitude  Plot  of  171  for  Unstable,  Nonminimum 
Phase  Syston  Example  (Small  Noise  Case) 


Frequency  (rad/sec) 


Figure  6*56.  1-DOF  Bode  Magnitude  Plot  of  171  for  Unstable,  Nonminimum 
Phase  System  Example  (Large  Noise  Case) 
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I¥Dqiiency  (rad/sec) 

Figure  6-57.  1-DOF  Bode  Magnitude  Plot  of  IGiH  for  Unstable, 

Nonminimum  Hiase  Syston  Example  (Small  Noise  Case) 


Frequency  (rad/Mc) 

Figure  6-58.  1-DOF  Bode  Magnitude  Plot  of  \GK\  for  Unstable, 

Nonminimum  Phase  System  Example  (Large  Noise  Case) 
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Figure  6*59.  l>DOF  Step  Re^xmses  for  Unstable,  Nonminimum  Phase 
Syston  Example  (Small  Noise  Case) 


Figure  6*60.  1-DOF  Step  Responses  for  Un^able,  Nonminimum  I^ase 
System  Example  (Large  Noise  Case) 
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6.4.2  2-DOF  Case.  The  plots  of  sensitivity  function  are  shovm  in  Figures 
6-61  and  6*62.  The  primary  difference  in  these  plots  from  the  other  2>DOF  examples 
is  that  die  overshoot  above  0  dB  is  much  higher,  indicating  worse  stability  margins. 

In  the  large  noise-small  disturbance  case,  the  DC  gain  of  |5|  is  about  13  dB,  showing 
how  small  |G/ir]|  must  be  at  low  frequency  to  accommodate  the  large  noise.  Figures  6- 
63  and  6-64  show  plots  of  the  complimentary  sensitivity  function.  In  die  large  noise- 
small  disturbance  case,  the  DC  gain  of  |7]  is  also  above  0  dB  (this  time  at  about  18 
dB).  The  other  trends  are  unchanged  from  the  previous  examples. 

Figures  6-6S  and  6-66  show  plots  of  jCTATil  and  \GK^  for  the  small  noise  case, 
and  Figures  6-67  and  6-68  show  plots  of  IGUTil  and  \GK'^  for  the  large  noise  case.  The 
only  significant  difference  between  these  plots  and  those  of  the  stable,  nonminimum 
phase  system  is  die  large  noise-small  disturbance  case  of  IGATjl.  In  this  case,  die  DC 
gain  of  the  function  is  at  about  -IS  dB  instead  of  the  normal  0  dB,  which  shows  the 
flexibility  of  the  shiqie  of  \GK^  in  taking  on  whatever  magnitude  is  necessary  to  track 
the  command  input.  In  both  the  large  and  small  noise  cases,  the  plots  fall  below  0  dB 
before  they  reach  1  rad/sec  so  that  no  energy  is  put  into  the  system  at  the  frequencies 
of  the  unstable  pole  or  zero. 

The  plots  of  \SGK‘^  and  the  2-DOF  step  response  are  identical  to  Figures  6-51 
and  6-52,  respectively. 

6.4.3  Comparison  Between  I-DOF  and  2~DOF  Cases.  The  unstable, 
nonminimum  phase  system  example  represents  the  most  challenging  of  the  SISO 
examples  presented.  Thus,  the  loop  shape  plots  exhibit  all  of  the  characteristics 
associated  with  right-half  plane  poles  and  zeros.  Both  the  1-DOF  and  the  2-DOF 
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systems  show  a  great  deal  of  overshoot  above  die  0  dB  line  in  the  sensitivity  and 
complimentary  sensitivity  plots.  As  in  the  other  types  of  examples,  when  the  assumed 
noise  is  large  compared  to  the  plant  disturbance  in  the  2-DOF  case,  die  DC  gain  of  S 
becomes  larger  (as  in  die  1-DOF  case),  and  in  some  cases,  it  is  above  0  dB  (which 
does  not  occur  in  the  1-DOF  case).  In  the  complimentary  sensitivity  plots,  the  2-DOF 
case  shows  that  as  the  assumed  noise  becomes  larger  relative  to  the  disturbance,  the 
rolloff  frequency  decreases  (as  in  the  1-DOF  case),  and  in  the  smallest  disturbance 
case,  the  DC  gain  of  |7]  is  above  0  dB  (\^ich  does  not  occur  in  the  1-DOF  case).  The 
most  noticeable  difference  between  the  1-DOF  and  2-DOF  system  is,  once  again,  in 
the  step  responses.  In  die  1-DOF  case,  the  responses  exhibit  increasing  initial 
undershoot  as  the  assumed  disturbance  level  increases  (due  to  the  nonminimum  phase 
zeros),  as  well  as  significant  overshoot  of  the  final  value  (due  to  the  unstable  poles). 
The  2-DOF  system,  however,  utilizes  the  ciqiabilities  of  the  prefilter  to  maintain 
tracking  over  a  wide  frequency  range  and  still  produces  the  optimal  step  response  for 
the  given  tracking  weight,  ^bich,  for  this  case,  exhibits  minimum  initial  undershoot 
and  no  overshoot  of  the  final  value. 
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Igurc  6<61.  2-DOF  Bode  Magnitude  Plot  of  LSI  for  Unstable,  Nonminimum 
niase  System  Example  (Small  Noise  Case) 


Figure  6>62.  2-DOF  Bode  Magnitude  Plot  of  LSI  for  Unstable,  Nonminimum 
Phase  System  Example  (Large  Noise  Case) 
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Figure  6*63.  2-DOF  Bode  Magnitude  Plot  of  171  for  Unstable,  Nonminimum 
Phase  System  Example  (Small  Noise  Case) 


Figure  6-64.  2-DOF  Bode  Magnitude  Plot  of  171  for  Unstable,  Nonminimum 
Phase  System  Example  (Large  Noise  Case) 
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Figure  6>65.  2-DOF  Bode  Magnibide  Plot  of  IGAT,!  for  Unstable, 

Nonminimum  I%ase  Systnn  Example  (Small  Noise  C^) 


Figure  6-66.  2-DOF  Bode  Magnitude  Plot  of  IGiiTjl  for  Unstable, 

Nonminimum  Phase  System  Example  (Small  Noise  Case) 
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Figure  6>67.  2-DOF  Bode  Magnitude  Plot  of  K/AT,!  for  Unstable, 

Nonminimum  I^ase  System  Example  (Large  Noise  Case) 


Figure  6*68.  2-DOF  Bode  Magnitude  Plot  of  IGATjl  for  Unstable, 
Nonminimum  I%ase  Example  (Large  Noise) 
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6.5  System  Type  Examples 

In  Ais  section,  the  effects  of  applying  different  inputs  to  both  1-DOF  and  2- 
DOF  ^sterns  are  investigated  to  find  out  udiedier  2-DOF  systems  offer  the  same 
advantages  diat  th^  do  vdien  die  input  is  only  a  step.  The  example  to  be  used  is  the 
unstable,  nonminimum  phase  plant  from  section  6.4,  looking  only  at  the  large  noise 
case.  This  plant  was  chosen  because  it  is  the  most  challenging  to  control,  and  the 
differences  between  the  1-DOF  and  2-DOF  time  responses  are  the  most  obvious 
(especially  in  the  large  noise  case). 

Recall  from  Chapters  II  and  m  that  system  type  is  defined  as  the  degree  of  the 
polynomial  input  for  which  the  steady  state  error  is  equal  to  a  non-zero  constant. 
Furthermore,  if  die  system  is  type  0,  then  it  produces  a  constant  error  to  a  step  input; 
if  the  system  is  type  I,  then  it  tracks  a  step  input  perfecdy  and  produces  a  constant 
error  to  a  ramp  input;  if  the  system  is  type  II,  then  it  tracks  step  and  ramp  inputs 
perfecdy  and  produces  a  constant  error  to  a  par^olic  input;  and  so  on.  These 
concepts  are  repeated  here  because  they  have  a  direct  bearing  on  vdiat  is  expected  in 
the  time  response  plots  for  the  aforementioned  example. 

Also,  recall  from  Chrqiter  V  that  the  tracking  weight,  in  die  H2 
optimization  setup  can  also  be  used  to  shrqie  the  tracking  behavior  of  the  system  with 
regard  to  ^stem  type.  If  resembles  an  integrator,  then  the  system  will  exhibit  type 
I  behavior;  if  it  resembles  a  ramp  input,  then  the  system  will  exhibit  type  n  behavior; 
etc.  The  shqiing  of  is  the  method  that  will  be  used  to  produce  the  various  system 
types  in  both  the  1-DOF  case  and  the  2-DOF  case.  In  both  cases,  the  number  of  low 
frequency  poles  near  the  origin  used  to  model  is  directly  proportional  to  the 
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number  of  low  frequracy  poles  near  die  origin  in  the  resulting  compensator.  It  should 
be  noted  that  in  die  2-DOF  case,  the  number  of  poles  at  or  near  the  origin  in  eidier 
or  ATj  cannot  necessarily  be  attributed  to  system  type  behavior,  in  general.  Since 
compensators  are  being  utilized,  however,  ifTj  and  K2  share  the  same  poles,  and  in  diis 
particular  case,  the  number  of  poles  at  or  near  die  origin  in  the  compensators  has  a 
direct  bearing  on  the  system  type  behavior.  For  comparison  purposes,  the  following 
values  for  the  tracking  weight  were  utilized  to  produce  die  various  system  type 
behaviors  (through  type  m): 


or  =  10 

^  s  +  0.001 

(Type  1) 

(6.5) 

Hr  ,  100 

^  (s  +  0.001)* 

{Type  11) 

(6.6) 

HT  =  1000 

{Type  III) 

(6.7) 

’  (s  *  0.001)* 

These  weights  are  plotted  in  Figure  6-69.  Since  the  only  performance  characteristic  of 
interest  regarding  system  type  is  the  ability  of  the  system  to  track  different  inputs,  the 
only  plots  that  will  be  presented  in  this  section  will  be  time  response  plots. 

6.5.1  Type  I  System.  The  examples  presented  in  this  chapter  in  sections  6.1 
through  6.4  were  essentially  type  I  systems.  In  all  cases,  the  tracking  weight  had  a 
pole  at  s  =  -0.01.  In  these  system  type  examples,  the  pole(s)  was  moved  a  bit  closer 
to  the  origin  (to  s  =  -O.OOi)  just  to  illustrate  the  effect  of  having  the  poles  of  the 
compensator  very  close  to  the  origin. 

Figures  6-70  and  6-71  show  the  step  responses  of  the  1-DOF  and  2-DOF 
systems,  respectively.  The  1-DOF  response  is  almost  identical  to  Figure  6-60,  in  that 
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Figure  6*69.  Bode  Magnitude  Plot  of  Tracking  Weights  for  Syston 
Type  Examples 


it  exhibits  extreme  nonminimum  phase  bdiavior  (initial  underdioot)  and  overshoot  of 
the  final  value.  The  trends  in  the  reqionse  for  the  increasing  disturbance  are  as 
previoudy  described  -  as  the  disturbance  size  increases,  the  initial  undershoot 
increases  and  die  settling  time  gets  faster.  Depmding  on  the  disturbance,  this  system 
settles  out  between  4  and  10  seconds,  with  its  oscillatory  motion  ranging  firom  -6  to 
44.25.  The  2-DOF  step  re^nse  exhibits  no  overshoot  of  the  final  value,  settles  out 
in  about  3  seconds,  and  diows  an  initial  underdioot  of  -0.5.  Clearly,  the  2-DOF 
^stem  produces  a  much  more  favorable  step  re^nse  in  this  case  (note  scale  changes 
to  emphasize  differences). 

The  errors  to  a  ramp  input  are  diown  in  Figures  6-72  and  6-73.  These  errors 


6-55 


are  e}q>ected  to  reach  a  constant  non-zero  value  for  bodi  systems  since  they  are  botfi 
type  I.  The  2-DOF  system  steadies  out  in  under  three  seconds  at  a  value  of  +1.1.  The 
1-DOF  ^stem  steadies  out  between  6  and  25  seconds,  at  values  ranging  from  -2  to 
-13.  The  overshoot  for  all  three  disturbance  sizes  is  about  1.8. 
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Figure  6-72.  1-DOF  Error  Reposes  to  a  Ramp  Input 
(Type  I  System  Example) 


Figure  6-73.  2-DOF  Error  Response  to  a  Ramp  Input 
(Type  I  System  Sample) 
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6.5.2  T)^  U  System.  For  tiiis  case,  another  pole  near  the  origin  has  been 
added  to  the  tracking  weight,  producing  an  additional  equivalent  pole  in  the 
compensator  of  bodt  systems.  The  1-DOF  and  2-DOF  step  responses  are  shown  in 
Figures  6-74  and  6-7S.  In  die  1-DOF  case,  the  responses  all  settle  out  in  about  4-S 
seconds,  i^ich  is  a  marked  improvement  over  the  type  I  case,  but  the  oscillatory 
motion  is  much  greater,  especially  in  the  large  disturbance  case.  It  ranges  between  -30 
and  +22.  The  2-DOF  response  shows  the  same  trend  by  showing  a  quicker  setding 
time  of  2.S  seconds  and  more  oscillatory  motion  than  in  the  type  I  case,  ranging  from 
4.1  to  +2.5. 

The  ramp  error  responses  are  shown  in  Figures  6-76  and  6-77.  Because  these 
are  now  type  n  systems,  diey  now  track  ramp  commands  perfectly.  Depending  on  the 
disturbance  level,  the  1-DOF  system  setdes  out  at  zero  between  5  and  15  seconds, 
with  oscillatory  motion  ranging  between  -2.1  and  +3.2.  The  2-DOF  ^stem  setdes  out 
at  zero  in  about  3  seconds,  with  an  initial  overshoot  of  +0.9. 

The  parabola  error  responses  are  shown  in  Figures  6-78  and  6-79.  As 
expected,  both  ^sterns  produce  constant  errors  to  a  parabolic  inputs.  The  1-DOF 
system  setdes  out  between  4  and  15  seconds  at  values  between  -16.5  and  -2.2.  The 
small  disturbance  case  is  significandy  slower  and  has  a  larger  error  than  the  two  larger 
disturbances.  The  2-DOF  system  setdes  out  in  about  3  seconds  at  a  value  of  1.3,  with 
no  initial  undershoot. 
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Figure  6*74.  1-DOF  Step  Refuses  for  Type  n  Syston  Example 


Figure  6-75.  2-DOF  Step  Re^onse  for  Type  n  System  Example 
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Figure  6>7<.  1-DOF  Error  Responses  to  a  Ramp  Riput 
(Type  n  Systnn  Example) 
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Figure  6*77.  2-DOF  Error  Response  to  a  Ramp  Riput 
(Type  n  System  Example) 
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Figure  6*78.  1-DOF  Error  Refuses  to  a  Parabolic  Input 
(Type  n  Syston  Example) 


Figure  6-79.  2-DOF  Error  Response  to  a  Parabolic  Input 
(Type  n  Systran  Example) 
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6,5,3  Type  III  System.  For  this  final  case,  there  are  now  three  poles  near  the 
origin  in  tiie  tracking  weight,  which  produce  diree  poles  near  the  origin  in  the  1-DOF 
and  2-DOF  compensators.  Both  systems  should  be  able  to  track  step,  ramp,  and 
parabolic  inputs  with  zero  steady  state  error.  The  1-DOF  and  2-DOF  step  responses 
are  shown  in  Figures  6-80  and  6-81.  The  settling  times  continue  to  get  faster,  but  die 
amount  of  overshoot  also  increases,  especially  in  the  1-DOF  case.  The  1-DOF  setding 
times  all  fall  under  two  seconds,  v^ile  the  overshoot  for  the  large  disturbance  is  ±130. 
The  2-DOF  setding  time  is  about  one  second,  with  oscillatory  motion  ranging  between 
-20  and  +16. 

The  ramp  error  responses  are  shown  in  Figures  6-82  and  6-83.  The  1-DOF 
system  setdes  out  between  4  and  10  seconds  (depending  on  the  disturbance  level), 
with  the  overshoot  at  ±12.  The  2-DOF  system  setdes  out  in  2.5  seconds  with 
oscillatory  motion  ranging  between  -1.1  and  +2.4. 

The  parabola  error  responses  are  shown  in  Figures  6-84  and  6-85.  The  two 
larger  disturbances  in  the  1-DOF  case  produced  a  setding  time  of  about  6  seconds,  but 
the  smallest  disturbance  case  resulted  in  the  ^stem  going  unstable,  due  to  numerical 
errors.  If  there  had  been  a  true  triple  pole  at  the  origin,  then  the  system  would  have 
produced  zero  steady  state  error.  However,  the  fV^  tracking  weight  only  contains  a 
close  approximation  to  a  triple  pole  at  the  origin,  and  the  truncation  error  caused  by 
the  approximation  eventually  affects  the  system  response,  especially  vsdien  the 
measurement  noise  is  large  relative  to  the  plant  disturbance.  A  better  approximation 
would  also  be  to  make  the  DC  gain  of  larger.  The  2-DOF  error  response  settles  out 
in  three  seconds,  with  an  initial  overshoot  of  +1.1. 
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It  is  interesting  to  note  that  for  the  small  disturbance  case  for  the  type  III 
system,  the  1-DOF  step  r^ponse  looks  a  lot  like  the  2-DOF  step  response  with  respect 
to  the  settling  time  and  amount  of  oscillatory  motion.  When  the  1-DOF  system  is 
given  ramp  and  parabolic  inputs,  however,  the  responses  become  much  worse.  The 
reason  for  this  goes  back  to  the  idea  that  the  tracking  weight  does  not  really  have  a 
triple  pole  at  the  origin.  From  Figure  6-69,  the  type  m  tracking  weight  does  have  a 
large  amount  of  energy  at  low  frequency,  enabling  a  favorable  1-DOF  step  response. 
As  previously  mentioned,  vdien  the  input  is  higher  order,  the  1-DOF  responses  rapidly 
deteriorate. 
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Figure  6*80.  1-DOF  Step  Responses  for  Type  m  Syatan  Example 


Figure  6*81.  2-DOF  Step  Req}onse  for  Type  m  System  Example 
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Figure  6>82.  l>DOF  Error  Reposes  to  a  Ramp  Input 
(Type  m  System  Example) 
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Figure  6-83.  2-DOF  Error  Response  to  a  Ramp  Input 
(Type  in  System  Example) 
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Figure  <>-84.  1-DOF  Error  Reposes  to  a  Parabolic  I^put 
(Type  in  Syston  Example) 


Figure  6-85.  2-DOF  Error  Response  to  a  Parabolic  Input 
(Type  m  System  Example) 
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6,6  Low  Frequency  Measurement  Noise  Example 


In  all  of  tfte  examples  presented  dius  far,  both  the  1-DOF  and  2-DOF  ^sterns 
have  beoi  modeled  with  plant  disturbances  that  contain  definite  low  frequency  content 
(Figure  5-3)  and  measurement  noises  that  contain  definite  high  frequency  content 
(Figure  5-4).  One  of  die  basic  ideas  emphasized  in  Chapter  V  is  diat  the  system  loop 
sh^e,  whether  1-DOF  or  2-DOF,  must  exhibit  certain  properties  at  low  frequencies  to 
deal  widi  commands  and  disturbances,  and  odier  properties  at  high  frequencies  to  deal 
with  measurement  noises.  The  idea  that  die  measurement  noise  m^  contain  frequmcy 
content  in  die  same  range  that  commands  and  disturbances  fall  has  not  been 
considered.  This  section  of  the  chapter  will  explore  die  effects  of  low  frequency 
measurement  noise  on  system  loop  shape  and  die  corresponding  time  responses. 

The  plant  to  be  used  for  this  example  is  once  again  the  unstable,  nonminimum 
phase  i^stem  from  section  6.4.  The  measurement  noise  weight  fi’om  the  H2 
optimization  setup,  W^,  must  be  modified  to  overlap  die  frequency  range  that  is  used 
to  model  die  plant  disturbances.  The  modified  measurement  noise  weight  is  given  by 

jy  _  20  (s  0.0001)  /g  gx 

’  ~  s  *  0.001 

The  Bode  magnitude  plot  of  this  weight  compared  to  the  medium  sized  disturbance 
weight  (widi  a  DC  gain  of  20  dB)  is  shown  in  Figure  6-86.  The  medium-sized 
disturbance  weight  is  shown  on  the  plot  only  to  show  the  overlap  of  the  frequencies. 
Ail  three  disturbance  weights  will  be  reflected  in  the  loop  shiqie  and  time  response 
plots  for  this  example.  The  frequency  and  time  response  plots  generated  from  this 
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Figure  6>86.  Bode  Magnitude  Hot  of  Low  Frequmcy  Measurraiait 
Noise  Weight 


example  are  intraded  to  be  compared  to  the  large  noise  plots  from  die  unstable, 
nonminimum  phase  example  from  section  6.4. 

For  die  1-DOF  case,  traditional  loop  dicing  theory  states  that  tCAH  must  be 
large  at  low  frequency  for  good  disturbance  rejection  and  small  at  high  frequency  for 
good  measuranent  noise  rejection.  The  previous  examples  have  shown  that  when  die 
measurement  noise  is  large  relative  to  the  disturbance,  the  size  of  \GK\  is  drivoi  down 
over  all  frequencies.  Thus,  in  the  small  disturbance  case,  \GK\  is  expected  to  be  small 
overall,  regardless  of  the  frequency  range  of  the  noise.  When  the  disturbance  is  larger, 
the  size  of  KSKl  is  not  as  predictable.  In  general,  it  should  be  expected  diat  die  size  of 
\GK\  will  be  drivoi  by  the  low  frequency  measurement  noise;  therefore,  \GK\  should  be 
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driven  down  at  the  frequency  where  the  noise  first  becomes  significant.  In  any  case, 
the  time  responses  should  be  worse  than  the  case  where  die  noise  is  restricted  to  die 
high  frequency  region. 

For  die  2-DOF  case,  there  is  still  a  potential  conflict  between  disturbance 
rejection  and  measurement  noise  rejection,  in  diat  they  both  depend  on  the  size  of 
IGiTil.  When  the  disturbance  is  small,  there  should  be  no  problem  for  the  2-DOF 
system  to  accomplish  bodi  objectives  because  does  not  need  to  be  large  to  reject 
the  disturbance,  and  it  should  be  small  to  reject  the  noise  anyway.  |GAr2l  compensates 
for  the  reduced  size  of  |G^i|  by  shaping  itself  in  such  a  way  that  tracking  is 
maintained  over  as  large  a  frequency  range  as  possible.  For  the  large  distinbance-large 
noise  case,  however,  there  is  still  a  conflict  between  die  disturbance  and  noise 
rejection  objectives.  From  a  logical  standpoint,  it  appears  that  diis  could  result  in  a 
worse  time  response  dian  the  case  vdiere  the  noise  contains  only  higher  frequencies. 

6.6.1  l~DOF  Results.  Figure  6-87  shows  die  1-DOF  sensitivity  function  for 
die  low  frequency  measurement  noise  case.  It  is  immediately  apparent  that  in  die  two 
smaller  disturbance  cases,  |5|  is  at  0  dB  between  0.0001  and  0.1  rad/sec,  which 
suggests  diat  \GK\  is  very  small  in  this  region.  Recalling  die  steady  state  error  barrier 
from  Figure  2-3,  it  does  not  ^pear  that  die  functions  in  the  two  smaller  disturbance 
cases  resemble  "good”  sensitivity  functions.  |5|  for  the  largest  disturbance  resembles  a 
more  conventional  sensitivity  function. 

Figure  6-88  shows  a  plot  of  the  complimentary  sensitivity  function.  Again,  the 
functions  for  die  two  smaller  disturbance  cases  exhibit  undesirable  behavior  in  that 
diey  "droop"  below  the  0  dB  line  at  frequencies  before  rolling  off,  indicating  that  |GA| 
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is  small  in  that  region.  Figure  6-89  verifies  die  fact  that  \GK\  is  indeed  small  in  the 
frequMicy  region  that  |71  exhibits  the  droop  (in  the  small  disturbance  cases).  Hie  DC 
gain  of  |GX]  is  also  very  low  for  those  two  cases  compared  to  the  large  disturbance 
case. 

The  step  responses  in  Figure  6-90  indicate  diat  when  the  assumed  disturbance 
is  large,  the  response  is  virtually  unchanged  from  the  case  where  the  measurement 
noise  is  only  at  high  frequency  (Figure  6-60).  When  the  disturbance  is  smaller, 
however,  shqie  of  the  response  does  not  change  much,  but  die  steady  state  value  does 
not  even  reach  unity.  The  reason  for  diis  is  diat  is  very  small  in  the  region  i^ere 
die  disturbance  and  measurement  noise  frequencies  overly,  which  is  where  it  needs  to 
be  large  to  produce  a  favorable  response.  It  should  be  emphasized  here  diat  in  die 
large  disturbance  case,  the  step  response  showed  perfect  tracking  performance  using 
die  evaluation  model  from  Chtqiter  V,  which  has  no  noises  acting  on  the  system. 

Since  the  assumed  large  disturbance  causes  \GK\  to  be  large  at  low  frequency,  and  |71 
is  at  0  dB  in  the  frequency  range  \^ere  the  measurement  noise  first  becomes  large, 
the  system  would  track  both  the  command  input  and  the  noise  in  that  frequency  range. 
In  fact,  if  the  evaluation  model  were  altered  so  that  low  frequency  measurement  noise 
acted  on  the  system,  dien  the  step  response  would  show  the  amplified  noise  in  it. 
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Figure  6-87.  1-DOF  Bode  Magnitode  Plot  of  ISI  for  Low  Frequency 
Measuremoit  Noise  Example 


Figure  6-88.  1-DOF  Bode  Magnitude  Plot  of  171  for  Low  Frequency 
Measurement  Noise  Example 
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Figure  6*89.  l-DOF  Bode  Magnitude  of  IGATt  for  Low  Frequency 
Measurement  Noise  Example 


Figure  6-90.  l-DOF  Step  Refuses  for  Low  Frequency  Measurement 
Noise  Rumple 


6-73 


6.6.2  2-DOF  Resuits.  Figures  6-91  and  6-92  show  die  2-DOF  sensitivity  and 
complimoitaiy  soisitivity  functions.  For  the  small  disturbance  cases,  |iS]  is  at  0  dB 
betweoi  0.0001  and  0.1  rad/sec,  and  |71  exhibits  a  "droop"  in  diis  frequency  range.  In 
die  large  disturbance  cases,  liSI  and  |7]  more  closely  resemble  re^ective  "good"  shiqies. 
Figures  6*93  and  6-94  show  the  magnitude  plots  of  \GKi\  and  \GK^.  In  the  large 
disturbance  case,  exhibits  high  gain  at  low  frequency  and  low  gain  at  high 
frequency  and  is  virtually  matched  by  [GK^^.  In  the  two  smaller  disturbance  cases, 
however,  IGilTil  has  the  same  droop  that  |7]  does,  but  \GK‘^  is  at  0  dB  in  this  frequency 
range.  The  plot  of  j-SG/TJ  in  Figure  6-95  shows  that  die  2-DOF  system  tracks 
commands  perfectiy  up  to  the  bandwidth  frequency  of  10  rad/sec.  The  step  reqxmse 
in  Figure  6-96  shows  that  even  with  a  significant  low  frequency  measurement  noise, 
die  2-DOF  system  still  produces  the  same  step  response  with  die  noiseless  evaluation 
model  as  die  case  vidiere  the  noise  contains  stricdy  high  frequencies  (Figure  6-S2). 

6.6.3  Comparison  Between  1-DOF  and  2-DOF  Cases.  The  1-DOF  and  2- 
DOF  sensitivity  and  complimentaiy  sensitivity  function  plots  are  virtually  identical. 

The  primary  difference  between  die  1-DOF  and  2-DOF  ^sterns  comes  in  the  shape  of 
the  loop  transfer  functions.  |G^  in  the  1-DOF  system  and  \GKi\  in  die  2-DOF  ^stem 
are  almost  identical  in  that  the  small  disturbance  cases  still  have  the  low  frequency 
droop.  In  die  2-DOF  case,  however,  jG^ji  remains  at  0  dB  throughout  die  range  of 
the  droop,  diereby  maintaining  its  command  following  bandwidth.  As  previously 
mentioned,  the  1-DOF  step  responses  in  die  small  disturbance  cases  are  un^le  to  track 
commands  widiout  producing  a  large  steady  state  error,  vdiile  the  2-DOF  system  does 
not  show  any  error.  The  1-DOF  step  response  in  die  large  disturbance  case  does  not 
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show  any  steady  state  error,  but  its  initial  undershoot  is  significantly  larger  than  diat  of 
die  2-DOF  case.  In  the  large  disturbance  case,  the  2-DOF  step  response  is  subject  to 
die  same  caveat  as  die  1-DOF  response  --  diat  is.  since  the  system  tracks  commands 
perfecdy  in  the  0.0001  to  0.1  rad/sec  frequency  range,  it  will  also  track  measurement 
noise  in  diat  frequency  range.  Thus,  i^en  the  assumed  disturbance  is  large,  the 
resulting  1-DOF  and  2-DOF  controllers  would  show  significant  amounts  of  noise  in 
their  step  responses  if  there  is  low  frequency  measurement  noise  in  the  evaluation 
model.  When  die  assumed  disturbance  is  smaller,  however,  neither  type  of  system 
allows  die  low  frequency  noise  to  pass  dirough  to  the  ouqiut.  The  penalty  for  the  1- 
DOF  system  in  this  case  is  that  it  is  also  unable  to  track  commands  and  shows  a  large 
steady  state  error  to  a  step  input.  The  2-DOF  system  compensates  for  the  reduced 
loop  gain  by  reshtqiing  the  K2  controller  to  maintain  perfect  tracking  and  shows  no 
steady  state  error  to  a  step  command. 
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Figure  6-91.  2-DOF  Bode  Magnitude  Plot  of  LSI  for  Low  Frequency 
Measuiemoit  Noise  Example 


Figure  6-92.  2-DOF  Bode  Magnitude  Plot  of  171  for  Low  Frequency 
Measuranent  Noise  Example 
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Figure  6*93.  2'DOF  Bode  Magnitude  Plot  of  IGiir,l  for  Low  Frequency 
Measurement  Noise  Example 


Figure  6-94.  2-DOF  Bode  Magnitude  Plot  of  IGiiTjl  for  Low  Frequency 
Measurement  Noise  Example 
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Figure  6*95.  2-DOF  Bode  Magnitude  Plot  of  l5GAr2l  Frequency 

Measurement  Noise  Example 


Figure  6-96.  2-IX)F  Step  Re^nse  for  Low  Frequency  Measurement 
Noise  Sample 
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6.y  Chapter  Summary 

Altfiough  this  chapter  makes  heavy  use  of  graphs  to  illustrate  tfie  similarities 
and  differences  between  1-DOF  and  2-DOF  controllers,  there  are  several  generalities 
that  can  be  extracted  from  the  information.  Widi  regard  to  stability  margins,  it  is 
obvious  from  the  graphs  of  |iS]  and  |7]  that  unstable  poles  cause  a  marked  decrease  in 
the  margins,  but  nonminimum  phase  zeros  have  an  even  more  detrimental  effect  on 
them.  The  most  important  conclusions  to  draw,  however,  deal  with  the  ability  of  die 
two  ^sterns  to  track  commands.  In  die  1-DOF  cases,  the  following  trends  were 
observed  (using  the  stable,  minimum  phase  ^stem  as  a  baseline): 

i)  In  all  cases,  the  rise  times  and  setding  times  improve  as  die  size  of  die  plant 
disturbance  is  increased.  These  two  times  become  slower  when  going  from  a  small 
measurement  noise  to  a  large  one. 

ii)  The  effect  of  unstable  poles  on  die  step  response  is  to  increase  the 
overshoot  of  the  final  value.  In  general,  the  amount  of  overshoot  decreases  as  the 
disturbance  increases  in  magnitude.  The  settling  times  are  also  slighdy  slower. 

iii)  The  primary  effect  of  nonminimum  phase  zeros  is,  of  course,  to  cause  the 
step  response  to  initially  head  in  a  negative  direction.  The  magnitude  of  this 
undershoot  increases  and  the  setding  time  gets  faster  as  the  disturbance  size  increases. 
A  large  measurement  noise  decreases  the  magnitude  of  the  undershoot  and  slows  die 
settling  time. 

iv)  In  the  unstable,  nonminimum  phase  system,  all  of  the  above  trends  are 
evident  in  the  step  responses,  showing  that  the  trends  are  not  mutually  exclusive.  In 
fact,  this  type  of  example  serves  as  a  good  model  to  illustrate  the  differences  between 
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the  I'DOF  and  2-  DOF  approaches  because  it  is  a  "worst  case”  scenario. 

In  all  of  the  classes  of  examples,  the  2-DOF  cases  only  produced  a  single  step 
response.  This  is  because  the  H2  optimal  2-DOF  controller  produces  the  step  response 
that  is  the  "best”  overall  for  a  given  tracking  weight,  no  matter  what  dre  size  of  the 
plant  disturbance  or  measurement  noise  is.  For  the  tracking  weight  given  in  Ch^ter 
V,  the  following  observations  were  made  regarding  the  2-DOF  step  responses:  (1)  In 
the  minimum  phase  case,  it  gave  the  fastest  response  that  does  not  exhibit  any 
overshoot;  (2)  In  the  nonminimum  phase  case,  it  gave  the  response  with  the  lowest- 
sized  initial  undershoot,  fastest  settling  time,  and  no  overshoot  of  the  final  value. 

With  all  of  diese  constraints,  the  2-DOF  response  may  not  always  be  the  most 
practical  solution  in  everyday  use.  For  instance,  if  restrictions  are  relaxed  to  allow  a 
10%  overshoot,  then  die  rise  time  could  become  quite  a  bit  faster.  If  die  2-DOF 
tracking  weight  is  adjusted  to  roll  off  at  a  higher  frequency,  then  die  resulting  step 
response  would  exhibit  a  faster  rise  time  widi  more  overshoot.  The  1-DOF  controllers 
also  allow  for  this  possibility,  but  the  problem  with  the  1-DOF  controllers  is  that  they 
are  susceptible  to  changes  in  disturbance  and  measurement  noise.  The  2-DOF 
approach  guarantees  a  response  widi  the  aforementioned  characteristics  regardless  of 
the  size  of  the  disturbance  or  measurement  noise. 

The  system  type  examples  show  that  even  with  different  types  of  inputs,  the  2- 
DOF  system  still  produces  better  tracking  results  over  the  spectrum  of  plant 
disturbance  and  measurement  noise  sizes.  The  1-DOF  system  is  prone  to  exhibit  more 
oscillatory  motion  with  poles  near  the  origin  in  the  compensator  than  the  2-DOF 
system  is.  The  2-DOF  system  still  exhibits  the  "best"  overall  responses  with  higher 
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order  inputs. 

The  low  frequency  noise  example  perhaps  best  shows  the  advantages  of  using 
2-DOF  controllers.  The  plots  indicate  diat  when  die  measurement  noise  contains 
energy  in  die  same  range  as  the  plant  disturbance,  the  noise  is  more  dominant,  and  die 
loop  gain  will  be  low  to  lessen  its  effects.  In  the  2-DOF  case,  the  prefilter 
compensator  makes  up  for  the  smaller  size  of  |G^||  by  shaping  itself  to  maintain  a 
high  bandwidth  frequency  for  command  following.  The  1-DOF  ^stem  simply  has  no 
mechanism  to  put  forth  high  gain  at  low  frequency,  yet  avoid  the  effects  of  low 
frequency  measurement  noise. 

The  bottom  line  to  this  chapter  is  that  2-DOF  controllers  exhibit  much  better 
tracking  characteristics  dian  1-DOF  controllers  for  a  variety  of  classes  of  plants  and 
system  types.  The  next  chiqiter  examines  the  differences  between  die  two  systems  ir.  a 
MEMO  example. 
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VIL  Tail^Controlled  Missile  Example  (MIMO) 


As  previously  mentioned  in  Chi4)ters  n  and  m,  loop  shaping  concepts  can  also 
be  applied  to  die  multivariable  case.  The  primary  change  is  diat  restrictions  on  die 
sizes  of  various  loop  transfer  functions  are  now  placed  on  the  minimum  or  maximum 
singular  values  of  the  function.  The  example  presented  in  this  chapter  is  an  autopilot 
for  a  tail-controlled  missile  and  is  taken  from  [1:16-17].  It  is  a  point  design  in  die 
flight  envelope  of  die  missile  and  is  not  intended  to  be  a  design  used  for  operational 
purposes.  The  purpose  of  diis  example  is  to  once  again  show  die  differences  in  loop 
shapes  and  time  responses  between  1-DOF  and  2-DOF  systems. 

The  linearized  model  of  the  missile  plant  consists  of  eight  states,  given  by 


5p 

fin  actuator  state  {roll) 

6^ 

fin  actuator  state  [pitch) 

8r 

fin  actuator  state  [yaw) 

a 

pitch -plane  angle  of  attack 

P 

aw -plane  angle  of  attack 

P 

roll  rate 

q 

pitch  rate 

r 

yaw  rate 

The  inputs  consist  of  roll-fin,  pitch-fm,  and  yaw-fin  deflections.  The  plant  state  space 
matrices  are  given  by 


7-1 


-4.0000£+02 
0.0000£+00 
0.0000£+00 
-8.7349£-01 
-2.461 1£ -01 
3.3304£+06 
-1.0950£+03 
-1.2044£+02 


0.0000£+00 

O.OOOOF+00 

0.0000£+00 

-9.2340£-02 

-9.8435£-01 

2.6737£+02 

-6.9672£-01 

1.9462£+02 


0.0000£+00 

-4.0000£+02 

0.0000£+00 

1.3044£+02 

6.2620£-01 

1.2009£+06 

1.2844£+05 

5.6387£+02 

0.0000£+00 
O.OOOOE +00 
0.0000£+00 
0.0000£+00 
0.0000£  +00 
-1.5949£+00 
0.0000£+00 
0.0000£  +00 


0.0000£  +00 
0.0000£+00 
-4.0000£  +02 
-4.5412£-01 
-1.2753£+02 
-4.2461£+03 
5.8025£+01 
1.2598£+05 

0.0000£+00 
O.OOOOE +00 
0.0000£+00 
1.0000£+00 
0.0000£+00 
0.0000£+00 
-1.5498£+00 
0.0000£+00 


0.0000£+00 
0.0000£+00 
0.0000£+00 
-9.8435£-01 
-9.2340£  -02 
-2.6737£+02 
-1.9462£+02 
6.9673£  -01 

0.0000£  +00 
O.OOOOE +00 
0.0000£+00 
0.0000£+00 
-1.0000£+00 
0.0000£+00 
0.0000£+00 
-1.5498£  +00 


1  0  0 
0  1  0 
0  0  1 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 


(7.3) 


0.0000£+00 

0.0000£+00 

0.0000£+00 

-5.4953£+02 

-1.9504£+03 


0.0000£+00 

0.0000£+00 

0.0000£+00 

1.3982£+03 

2.9126£+05 


0.0000£+00 

0.0000£+00 

0.0000£+00 

-2.8476£+05 

-1.0140£+03 


0.0000£+00 

0.0000£+00 

0.0000£+00 

-2.0619£+02 

-2.1979£+03 


(7.4) 


0.0000£+00 
0.0000£+00 
0.0000£+00 
-2.1979£+03 
-2.061 9£ +02 


1.0000£+00 
0.0000£+00 
0.0000£+00 
O.OOOOE +00 
0.0000£+00 


0.0000£+00  0.0000£+00 
I.0000£+00  0.0000£+00 
0.0000£+00  1.0000£+00 
0.0000£+00  0.0000£+00 
0.0000£+00  0.0000£+00 
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The  ^stem  has  no  transmission  zeros,  and  the  poles  are  givoi  by 

-1.5949 

-1.2209  +  13.9227 
-1.2209  -  13.9227 

-1.3133  +  13.9747  (7.5) 

-1.3133  -  13.9747 
-400 
-400 
-400 

The  plant  outputs  available  for  measurement  (from  the  C  matrix)  are  die  angular  rate 
in  each  of  the  three  axes  (p,  q,  and  r)  and  the  acceleration  in  each  of  the  two  lateral 
axes  {Ny  and  N^.  The  only  outputs  of  interest  in  this  example  are  the  two 
accelerations  and  the  roll  rate.  Specifically,  the  control  objective  is  to  have  the  system 
track  a  10  g  A/j  command,  widi  Ny  and  p  commanded  to  zero. 

Since  diis  is  a  MEMO  example,  the  loop  shiqpes  are  now  represented  by 
singular  value  plots  instead  of  single-line  Bode  magnitude  plots.  Thus,  it  would  be 
difficult  to  observe  trends  with  multiple  singular  value  plots  on  a  single  gr^h.  To 
alleviate  diis  situation,  only  plots  resulting  from  die  smallest  and  largest  plant 
disturbances  will  be  presented  in  this  chiqiter.  Furdiermore,  only  the  minimum  and 
maximum  singular  values  for  each  disturbance  will  be  plotted.  The  singular  values  for 
each  disturbance  will  be  plotted  with  a  single  line  type,  and  the  disturbances  will  still 
be  identified  by  dieir  DC  gain.  In  this  manner,  two  singular  value  plots  can  be  plotted 
together  on  one  griq}h,  and  the  trends  can  still  be  easily  observed.  As  in  the  odier 
examples,  die  plots  in  which  small  measurement  noises  were  used  will  be  on  the  top 
of  the  page,  and  the  ones  in  which  large  noises  were  used  will  be  on  the  bottom. 
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7.1  1~D0F  Setup  and  Results 

The  first  problem  immediately  encountered  when  solving  diis  example  is  the 
fact  that  diere  are  five  outputs  available  for  measurement,  but  only  three  commanded 
inputs.  Since  unity  feedback  dieory  requires  diat  die  number  of  commanded  inputs 
equal  the  number  of  ouqiuts  being  fed  back,  then  the  problem  arises  as  to  what  to  do 
with  die  remaining  two  plant  outputs.  A  S  x  3  selector  matrix  can  be  placed  in 
judicious  locations  in  the  1-DOF  optimization  block  diagram  from  Figure  S-1  in 
order  to  make  maximum  use  of  die  outputs  in  the  controller  design.  As  it  turns  out, 
after  several  iterations  of  placing  the  selector  matrices  in  different  locations,  die  two 
remaining  ouqiuts  {q  and  r)  have  no  bearing  on  the  system’s  tracking  performance.  It 
is  generally  not  a  valid  assumption  to  ignore  plant  outputs  that  are  not  commanded, 
but  that  assumption  is  valid  in  diis  particular  case.  Therefore,  for  the  purposes  of  the 
1-DOF  controller  designs  in  diis  chapter,  the  plant  C  matrix  will  be  truncated  to  reflect 
only  the  diree  outputs  of  interest.  This  setup  also  results  in  die  compensator  widi  the 
lowest  order,  which  is  desirable  from  a  practical  standpoint.  With  the  C  matrix 
modified  in  diis  manner,  the  system  has  transmission  zeros  located  at  s  -  +28.6,  s  = 
+2S.S,  s  =  -30.1,  and  s  =  -27.1,  indicating  that  this  is  a  nonminimum  phase  ^stem. 
Figure  S-1  is  still  valid  for  use  in  the  H2  algorithm,  die  only  difference  being  that  all 
of  the  signals  are  now  three-dimensional  vector  signals,  and  the  weight  blocks 
represent  3x3  transfer  fimction  matrices.  The  weights  are  the  same  as  those 
mentioned  in  Chapter  V,  except  they  now  have  been  multiplied  by  die  identity  matrix. 
Equations  (S.l)  and  (S.2)  are  also  still  valid  representations  of  the  input-to-output 
transfer  function  matrix  and  the  system  state  space,  respectively. 
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Figures  7-1  and  7-2  show  plots  of  die  MIMO  sensitivity  function,  |iS],  for  the 
small  and  large  noise  cases,  respectively.  Using  the  ideas  from  Ch^iter  U,  it  is 

eiqiected  diat  the  maximum  singular  value  of  S,  o[5],  should  be  small  at  low 
frequency.  As  the  disturbance  gets  larger,  o[5]  gets  smaller.  In  the  large  noise  case, 

these  same  trends  are  present,  except  in  the  small  disturbance  case,  o  [5]  is  not  as 
small  at  low  frequmcy  as  it  is  in  the  small  noise  case.  In  the  large  noise-large 
disturbance  case,  a[<S]  is  about  the  same  magnitude  as  it  is  in  the  small  noise-large 
disturbance  case,  but  the  minimum  singular  value  of  S,  o[>S],  is  much  lower  at  low 

frequency  in  die  small  noise-large  disturbance.  This  is  because  with  a  smaller  noise, 
the  system  could  potentially  exhibit  higher  loop  gain  at  low  frequoicy,  as  in  the  SISO 
case.  The  "hump"  above  0  dB  is  larger  in  the  small  noise  cases  dian  it  is  in  the  large 
noise  cases  (indicating  worse  stability  margins)  because  when  more  gain  is  put  into  a 
nonminimum  phase  system,  die  response  gets  faster,  but  the  ^stem  comes  closer  to 
instability. 

The  system  complimentary  sensitivity  function,  |71,  is  shown  in  Figures  7-3  and 
7-4  for  each  noise  size.  For  good  noise  rejection,  a  [7]  must  be  small  at  high 

frequency,  and  the  MIMO  bandwidth  is  defined  as  the  frequency  where  a  [7]  rolls  off. 

As  the  disturbance  size  increases,  the  bandwidth  frequency  gets  larger,  as  expected.  In 
the  large  noise  case,  the  trends  are  the  same,  but  the  bandwidths  are  significantly 
lower,  as  in  die  SISO  case.  Once  again,  in  the  small  noise  case,  there  is  a  larger 
increase  in  the  function  magnitude  before  rolling  off  than  in  the  large  noise  case 
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because  die  larger  loop  gain  causes  worse  stability  margins  in  diis  system. 

Figures  7-5  and  7-6  show  plots  of  die  loop  transfer  fimction  \GK\.  As 
described  in  Chapter  n,  o[(/^  must  be  large  at  low  {requmcy  for  good  tracking  and 

disturbance  rejection,  and  a[GK]  must  be  small  at  high  frequency  for  good 
measurement  noise  rejection.  These  trends  are  observed  in  the  figures.  As  die 
disturbance  size  increases,  a[GK\  also  gets  bigger,  as  does  the  separation  betwe«i 

g[GA^  and  o[Giir]-  The  trends  are  identical  in  both  die  small  and  large  noise  cases, 

but  in  die  large  noise  case,  the  the  plots  are  shifted  downward,  reflecting  die  reduced 
loop  gain. 

Figures  7-7  and  7-8  show  the  responses  to  a  10  g  step.  In  the  small  noise 
case,  die  response  setdes  out  slighdy  faster  in  die  large  disturbance  case  (O.S  seconds 
vs.  0.6  seconds),  but  diere  is  also  higher  overshoot  (bodi  positive  and  negative)  than 
die  small  disturbance  case.  In  the  large  noise  case,  die  large  disturbance  response  still 
setdes  out  in  less  than  a  second  (0.7  seconds),  but  the  small  disturbance  response  has  a 
setding  time  of  over  10  seconds.  The  overshoot  in  die  large  disturbance  case  is  much 
higher,  however.  In  both  noise  size  cases,  die  effect  of  increasing  die  disturbance  is  to 
decrease  the  setding  time,  but  at  the  same  time  increase  the  magnitude  of  the 
oscillatory  motion.  These  exact  trends  are  true  of  the  Ny  responses  (Figures  7-9  and  7- 
10)  and  the  roll  rate  responses  (Figures  7-11  and  7-12). 


Singular  Values  (dB)  Siqgnlar  Vakies  (dB) 


Figure  7'1. 


I'DOF  Singular  Value  Plot  of  LSI  for  Tail-Controlled  Missile 
Example  (Small  Noise  Case) 


Figure  7-2.  1-DOF  Singular  Value  Plot  of  ISI  for  Tail-Controlled  Missile 
Example  (Large  Noise  Case) 
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Singular  Vahus(dB)  Siqgolar  Values  (dB) 


10 


Figure  7-3. 


1-DOF  Singular  Value  Plcrt  of  171  for  Tail-ControUed  Missile 
Exan^ile  (Small  Noise  Case) 


Figure  7-4.  1-DOF  Singular  Value  Plot  of  171  for  Tail-Controlled  Missile 
Example  (Large  Noise  Case) 
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SiQgularValnes(dB)  ^  Sii^iilarV«liiet(dB} 


igure  7-5.  1-DOF  Singular  Value  Plot  of  }GK\  for  Tail-Cmitrolled  Misaile 
Example  (Small  Noise  Case) 


Frequency  (lud/^) 


Figure  7-6.  1-DOF  Singular  Value  Plot  of  HJK\  for  Tail-Controlled  Missile 
Example  (Large  Noise  Case) 
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Nonnal  Aocdention  (g‘8)  _  Nomul  Accdentioo  (g‘t) 


figure  7*7.  1>D0F  N,  Repose  to  a  10  g  Step  for  Tail-Controlled  Missile 
Example  (Small  Noise  Case) 


Figure  7-8.  1-DOF  N,  Re^nse  to  a  10  g  Step  for  Tail-Controlled  Misale 
Example  (Large  Noise  Case) 
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LatenlAoodeiati<n(g‘8)  UttnlAocdentioa(g‘8) 


Figure  7-9.  1-DOF  Nj  Req)onse  to  a  10  g  Step  in  N,  for  Tail-Controlled 
Missile  Exmple  (Small  Noise  Case) 


Time  (sec) 


Figure  7-10.  1-DOF  Re^nse  to  a  10  g  Step  in  N,  for  Tail-Controlled 
Missile  Example  (Large  Noise  Case) 
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Roll  Rate  (degree)  Rdl  Rite  (d^^) 


Figure  7«11.  1>D0F  Roll  Rate  Repose  to  a  10  g  Step  in  N,  for 

Tail-Controlled  Mi^le  Example  (Small  Noiae  Case) 


Figure  7-12.  1-DOF  RoU  Rate  Response  to  a  10  g  Step  in  N,  for 

Tail-Controlled  Missile  Example  (Large  Noise  Case) 
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7.2  2-DOF  Setup  and  Results 

Section  7.1  revealed  that  the  1-DOF  block  diagram  setup  has  some 
limitations  as  far  as  being  able  to  utilize  all  of  the  possible  plant  outputs  in  the  design 
of  the  compmsator.  The  2*DOF  setup  contains  a  bit  more  flexibility  than  the  1-DOF 
setup  does  because  the  K^  compensator  is  located  after  the  plant  in  die  control  loop 
instead  of  before  it.  In  this  manner,  all  five  of  die  plant  outputs  can  be  fed  direcdy 
into  instead  of  "throwing  away"  two  states.  The  Ki  compensator  then  becomes  a 
non-square  3  x  S  transfer  function  matrix.  The  modified  2-DOF  setup  is  shown  in 
Figure  7-13.  The  M  block  in  the  figure  is  die  selector  matrix  that  is  used  to  convert 
die  plant  output  signal  from  five  outputs  to  diree  and  is  given  by 


Af 


1  0  0  0  0 
0  0  0  1  0 
0  0  0  0  1 


(7.6) 


Figure  7-13.  2-DOF  MIMO  Optimization  Block  Diagram 
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The  numbers  located  at  various  points  around  the  figure  represent  the  number  of  inputs 
or  outputs  that  are  contained  in  die  respective  vector  signal.  Equation  (5.4)  is  still  a 
valid  representation  of  die  input-to-output  transfer  function  matrix,  except  for  the  last 
row,  the  Zj  output,  vdiich  needs  to  reflect  the  addition  of  die  selector  matrix.  The 
modified  transfer  function  matrix  is  given  by 


W^SGK^ 

W^[I  -  K^SG\K^ 

w^[i  -msgk;\ 


-fF^KjSfV^ 

-W^MSW^ 


-W,TW^  1 

r 

1  s 

'•c 

-W^K^SW^ 

W^MTW^ 

w 

(7.7) 


The  state  space  has  only  a  few  minor  modifications  from  equation  (5.6)  and  is  not 
shown  here. 

The  2-DOF  MEMO  sensitivity  function,  15],  is  shown  in  Figures  7-14  and  7-15. 
Since  the  quantity  GKi  is  die  product  of  a  5  x  3  matrix  and  a  3  x  5  matrix,  die 
resulting  matrix  is  5  x  5,  but  its  rank  is  only  three.  Therefore,  there  are  two  singular 
values  that  are  equal  to  zero  over  all  frequencies,  although  they  are  not  shown  in  the 

figures.  In  die  small  noise  case,  a[5]  is  very  similar  to  die  1-DOF  case,  in  that  it 
decreases  as  the  disturbance  size  increases.  In  the  large  disturbance  case,  however, 
g[5]  is  much  lower  than  it  is  in  the  1-DOF  case,  indicating  the  potential  for  higher 

loop  gain.  In  the  large  noise  case,  some  of  the  same  trends  are  present  as  in  the  SISO 
case,  such  as  the  singular  values  all  being  zero  when  the  disturbance  is  small,  showing 
the  Ki  compensator  minimizing  the  effects  of  the  noise. 

Figures  7-16  and  7-17  show  plots  of  die  MIMO  complimentary  sensitivity 
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function.  As  in  the  SISO  case,  the  bandwidth  increases  as  the  disturbance  size 
increases.  When  the  disturbance  is  large,  in  both  the  small  and  large  noise  cases, 

o[71  rises  above  0  dB  before  rolling  off,  indicating  deteriorating  stability  margins  due 

to  increased  loop  gain  in  a  nonminimum  phase  system.  In  die  small  noise>large 
disturbance  case,  die  singular  values  do  not  roll  off  until  after  10,000  rad/sec, 
compared  to  100  rad/sec  in  the  1-DOF  case.  In  the  large  noise-small  disturbance  case, 

o[7]  has  a  DC  gain  of  nearly  -100  dB,  again  showing  minimizing  the  effects  of 

the  noise. 

Figures  7-18  and  7-19  show  o[GA^i]  and  a[GK^  for  die  small  noise  case.  In 
the  large  disturbance  case,  the  separation  between  o[GA^,]  and  o[G/r,]  is  much  bigger 

dian  it  is  in  the  small  disturbance  case,  i^^ich  shows  the  potential  for  higher  loop  gain, 
as  expected  from  the  sensitivity  plots.  Also,  since  the  noise  is  so  small,  the  K2 
compensator  does  not  really  come  in  to  pliQr,  and  the  singular  values  of  GK^  and  GK2 
are  almost  die  same  over  all  frequencies. 

Figures  7-20  and  7-21  show  o[G/ir,]  and  cs[GK^  for  the  large  noise  case.  In 
the  large  disturbance  case,  the  trends  are  the  same  as  they  are  in  the  small  noise  case. 
In  the  small  disturbance  case,  however,  a[G^,]  is  held  down  by  the  dominant 

measurement  noise,  and  all  of  the  singular  values  of  GK2  are  at  0  dB  to  maintain 
tracking  performance,  as  in  the  SISO  case. 

Figure  7-22  shows  a  plot  of  ISG/Tjl,  which  is  the  same  for  all  of  the  plant 
disturbance  and  measurement  noise  combinations,  showing  the  bandwidth  frequency 
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for  tracking  commands.  Figure  7-23  shows  a  plot  of  the  step  response.  It  has  an 
initial  undershoot  of  less  dian  -1  and  a  settling  time  of  about  0.6  seconds,  with  no 
overshoot  of  die  final  value.  This  response  is  virtually  identical  to  die  small  noise- 
small  disturbance  case  for  the  1-DOF  system  and  represents  die  best  response  possible 
without  any  overshoot  of  the  final  value.  The  responses  of  Ny  and  roll  rate  (Figures  7- 
24  and  7-25,  respectively)  exhibit  responses  with  the  same  characteristics  of 
minimizing  die  over-  and  undershoot. 

Overall,  this  MIMO  example  served  to  reiterate  die  concepts  derived  in  Chapter 
m  and  those  learned  from  die  SISO  examples  in  Chapter  VI.  Whenever  the 
measurement  noise  is  large  enough  to  dominate  the  system,  the  2-DOF  system  is  able 
to  react  by  driving  down  the  singular  values  of  GK^  to  minimize  its  effects  and 
readjusting  die  singular  values  of  GK^  to  maintain  favorable  tracking.  The  1-DOF 
system  responses  are  particularly  susceptible  to  die  effects  of  large  measurement  noise, 
i^ich  serves  to  increase  the  overshoot  and  setding  time.  For  the  given  tracking 
weight,  die  2-DOF  system  still  produces  the  best  overall  time  responses  from  the 
standpoint  of  exhibiting  no  overshoot  of  the  final  value  and  minimizing  the  initial 
undershoot,  regardless  of  the  disturbance  or  measurement  noise  sizes.  For  this  MIMO 
example,  die  2-DOF  system  also  proves  to  make  use  of  all  five  available 
measurements,  vdiich  the  1-DOF  system  does  not.  In  a  more  complicated  MIMO 
example,  where  knowledge  of  certain  measurements  may  be  crucial  to  controlling 
other  commands,  the  2-DOF  system  would  be  at  a  distinct  advantage. 
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Singular  VahieaCdB)  Singular  Vahiea(dB) 


figure  7<14. 


2-DOF  Singular  Value  Plot  of  ISI  for  Tail-Controlled  Missile 
Example  (Small  Noise  Que) 


Figure  7-15.  2-DOF  Singular  Value  Plot  of  151  for  Tail-Controlled  Missile 
Example  (Large  Noise  Case) 
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•'igure  7>16.  2-IX)F  Singular  Value  Plot  of  171  for  Tail-Controlled  Missile 
Example  (Small  Noise  Case) 


Figure  7-17.  2-DOF  Singular  Value  Plot  of  171  for  Tail-Controlled  Missile 
Example  (Large  Noise  Case) 
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Figure  7«18.  2>D0F  Singular  Value  Plot  of  I6fi^,l  for  Tail-Controlled  Missile 
Example  (Small  Noise  Case) 


Figure  7-19.  2-DOF  Singular  Value  Plot  of  ICATjl  for  Tail-Controlled  Missile 
Example  (Small  Noise  Case) 
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gure  7«20.  2-DOF  Singular  Value  Plot  of  for  Tail-Controlled  Miaaile 
Exan^ile  (Large  Noiae  Caae) 


Figure  7-21.  2-DOF  Singular  Value  Plot  of  IGATjl  for  Tail-Ck)ntrolled  Missile 
Example  (Large  Noise  Case) 
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Figure  7<22.  2>DOF  Singular  Value  Plot  of  ISGAT,!  Tail-Qxitrolled 
Missile  Example  (any  size  noise) 


Figure  7-23.  2-DOF  N,  Re^nse  to  a  10  g  Step  for  Tail-Controlled  Missile 
Example  (any  size  noise) 
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Figure  7>24.  2-IX)F  N,  Repose  to  a  10  g  Step  in  N,  for  Tail-Controlled 
Missile  Ex^ple  (any  size  noise) 


Figure  7-25.  2-DOF  Roll  Rate  Re^nse  to  a  10  g  Step  in  N,  for 
Tail-Controlled  Missile  Example  (any  aze  noise) 
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VIIL  Conclusions  and  Recommendations 


8.1  Summary  attd  Concbtsions 

At  the  beginning  of  this  thesis,  it  was  stated  that  2-DOF  systems  were  assumed 
to  exhibit  better  tracking  qualities  than  1-DOF  systems,  and  two  of  the  main  qu^ons 
to  answer  were  how  much  and  under  wdiat  conditions  they  were  better.  The  approach 
to  quantifying  die  differences  between  die  two  systems  started  out  with  a  summary  of 
die  rules  governing  loop  shaping  in  the  1>D0F  case,  and  dirough  the  extension  of 
these  concepts  to  the  2-DOF  case,  key  differences  were  identified.  Several  types  of 
examples  were  then  completed  in  order  to  try  to  achieve  the  basic  control  objectives  of 
tracking  commands  and  rejecting  plant  disturbances  and  measurement  noises  utilizing 
both  types  of  systems. 

In  the  four  basic  SISO  examples,  the  frequency  re^nses  gave  the  basic  closed 
loop  shqies  that  could  be  expected  to  be  encountered  in  control  system  design. 
Although  the  nominal  plant  model  was  low  order  and  could  be  considered  very  simple, 
it  did  allow  for  the  effects  of  putting  a  pole  or  zero  into  die  right-half  plane  to  be 
easily  observed.  It  turned  out  that  nonminimum  phase  zeros  had  the  most  detrimental 
effect  on  the  tracking  ability  of  both  systems,  sharply  reducing  the  stability  margins 
and  causing  die  step  responses  to  initially  head  in  the  opposite  direction.  The  fact  that 
zeros  cannot  be  moved  through  feedback  control  is  a  universal  problem  in  control 
dieory.  Unstable  poles  also  had  adverse  effects  on  the  systems,  particularly  by 
reducing  the  stability  margins. 
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Most  importantly,  the  four  SISO  examples  revealed  die  conditions  under  which 
2-DOF  controllers  really  have  an  advantage  over  1-DOF  controllers.  In  cases  vdiere 
die  plant  disturbance  was  large,  the  time  responses  of  both  ^sterns  were  very  similar. 
In  fact,  die  1-DOF  responses  often  showed  faster  rise  and  setding  times.  The  one 
variable  which  brought  out  the  differences  was  measurement  noise.  Since  the  system 
reference  command  and  measurement  noise  are  related  to  the  output  by  the  same 
transfer  function  in  the  1-DOF  system,  it  is  readily  apparent  that  diere  is  always  a 
potential  conflict  between  command  following  and  noise  rejection.  The  noise  level 
could  be  viewed  as  a  "knob”  diat  cat  be  direcdy  tuned  to  degrade  the  tracking 
performance  of  the  system.  In  all  of  the  examples,  increased  measurement  noise 
caused  more  overshoot  and  longer  setding  times  in  the  1-DOF  step  responses.  The  2- 
DOF  controller's  structure  effectively  decouples  the  requirements  necessary  for  noise 
rejection  and  command  following  performance  by  allowing  the  feedback  controller  to 
concentrate  on  attenuating  measurement  noise  and  the  prefilter  controller  to  adjust 
accordingly  to  maintain  tracking  over  as  wide  a  frequency  range  as  possible. 

The  effects  of  measurement  noise  on  tracking  performance  was  made  more 
clear  in  the  low  frequency  measurement  noise  example.  When  the  measurement  noise 
frequencies  dip  into  the  region  normally  reserved  for  plant  disturbances  and  command 
inputs,  the  classic  tradeoff  between  tracking  and  noise  rejection  in  1-DOF  systems 
begins  to  break  down,  as  the  system  must  try  to  simultaneously  attenuate  the  noise 
while  at  the  same  time  place  energy  into  tracking  commands  -  an  impossible  task 
The  1  -DOF  step  response  for  this  example  showed  that  as  the  noise  level  increased,  so 
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did  the  steady  state  error  in  the  response,  due  to  the  fact  that  no  loop  gain  could  be 
applied  at  the  frequencies  v^ere  high  gain  is  normally  required  for  a  low  steady  state 
error.  The  2-DOF  system  alleviated  this  problem  by  shying  the  feedback  controller  to 
resemble  that  of  die  1-DOF  system,  thereby  holding  down  the  effects  of  the  noise. 
Meanwhile,  die  prefilter  controller  "picked  up  the  slack"  by  shying  itself  in  such  a 
way  as  to  preserve  its  tracking  characteristics.  In  fact,  the  2-DOF  system  acted  as  if  it 
could  not  tell  the  difference  between  high  or  low  frequency  measurement  noise. 
However,  if  the  evaluation  model  were  changed  so  that  low  frequency  measurement 
noise  actually  entered  the  system,  then  both  the  1-DOF  and  2-DOF  systems  would 
exhibit  "noisy"  responses  in  the  case  v^.ere  die  controllers  were  designed  with  both  a 
large  plant  disturbance  and  large  measurement  noise  in  mind. 

The  SISO  system  type  examples  demonstrated  one  of  the  defining 
characteristics  of  2-DOF  controller  tracking  performance.  That  is,  in  each  of  the 
examples,  the  2-DOF  system  produced  the  response  that  had  the  fastest  settling  time, 
subject  to  the  constraint  of  having  no  overshoot.  In  nonminimum  phase  systems,  the 
amount  of  initial  undershoot  was  also  minimi;Kd.  Even  the  error  responses  for  ramp 
and  parabolic  inputs  showed  the  2-DOF  controller's  response  to  have  very  little 
oscillatory  motion.  A  possible  criticism  of  the  limits  of  the  2-DOF  time  response  for 
die  given  tracking  weight  is  that  I -DOF  controllers  produce  responses  which  have 
faster  rise  and  settling  times  under  certain  conditions,  such  as  when  the  measurement 
noise  level  is  low  compared  to  the  disturbance  level.  That  perceived  drawback  can  be 
dismissed  when  the  fact  that  the  2-DOF  controller  produces  responses  with  the 
aforementioned  constraints  regardless  of  the  size  of  the  plant  distuioance  or 
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measurement  noise  is  considered.  The  tracking  weight  can  be  altered  to  contain  more 
high  frequency  contrat  in  order  to  improve  the  speed  of  the  transient  response  (and 
allow  more  overshoot).  In  addition,  it  is  possible  to  move  the  poles  of  the  prefilter  to 
different  locations  after  die  optimization  algorithm  is  complete,  thereby  improving  the 
response  without  affecting  the  disturbance  or  noise  rejection  properties  of  the  system. 

The  tail-controlled  missile  MIMO  example  reiterated  all  of  the  concepts  just 
discussed  for  die  SISO  case  and  brought  out  another  area  where  2-DOF  controllers 
may  offer  possible  advantages.  When  a  plant  is  modeled  as  a  non-square  system 
(different  number  of  outputs  than  inputs),  as  the  missile  in  this  example  was,  1-DOF 
systems  are  limited  insofar  as  how  they  make  use  of  all  available  information.  Since 
die  input  to  the  compensator  is  die  difference  between  the  command  inputs  and  the 
measured  outputs,  dien  the  number  of  outputs  used  by  the  optimization  algorithm  is 
limited  by  the  number  of  command  inputs.  In  2-DOF  systems,  on  the  other  hand,  the 
feedback  controller  is  located  after  the  plant  in  the  control  loop.  Therefore,  the 
feedback  controller  can  be  non-square  and  dius  make  use  of  all  of  the  plant  outputs  yet 
at  die  same  time  feed  back  only  the  i^ipropriate  number  of  outputs  to  be  summed  widi 
die  command  inputs. 

The  concept  of  a  2-DOF  controller  has  been  around  for  several  years,  but  it  is 
still  relatively  new  from  the  point  of  view  that  there  are  virtually  no  textbooks  that 
include  it  as  part  of  a  basic  discussion  of  linear  control  theory.  The  primary  purpose 
of  this  work  was  to  take  the  concept  of  a  2-DOF  controller  and  compare  it  to  the  basic 
1-DOF  model  as  far  as  vidiat  it  takes  to  achieve  control  objectives.  Through  the 
examples  presented,  it  is  hoped  that  it  has  become  clear  that  a  2-DOF  controller  offers 
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significant  design  advantages  over  the  conventional  1-DOF  model  over  a  wide  variety 
of  conditions. 


8,2  Recommendations  for  Future  Research 

The  2-DOF  controllers  employed  in  diesis  made  use  of  the  H2  optimization 
method,  which  has  the  limitation  of  constraining  the  feedback  compensator  and  the 
prefilter  to  share  the  same  poles.  There  has  been  some  work  done  on  designing  each 
compensator  separately  in  an  effort  to  optimize  each  one  without  being  limited  by  the 
other.  One  potential  research  topic  would  be  to  investigate  the  techniques  of  designing 
die  two  controllers  separately  and  comparing  diem  to  designs  vdiere  the  controllers  are 
designed  in  a  single  step.  Also,  since  the  2-DOF  system  is  a  basic  control  system 
model,  other  optimization  techniques  could  be  applied  to  it,  such  as  LQG,  and 
mixed  H2/HS,.  Finally,  it  would  be  interesting  to  see  just  what  DOF  ^sterns  are 
capable  of  in  practical  qiplications.  If  simulations  indicate  that  they  produce  superior 
time  responses  regardless  of  any  plant  disturbances  or  measurement  noises,  then  2- 
DOF  controllers  would  make  excellent  candidates  for  precision  control  applications. 
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